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OPTICAL WAVEGUIDE STRUCTURES 

CROSS-REFERENCE TO RELATED APPLICATIONS 

This is a Continuation-in-Part of United States patent application serial No. 
5 09/742,422 which itself was a Continuation-in-Part of United States patent application 
serial No. 09/629,816 filed July 31, 2000, now United States patent No. 6,442,321 
issued August 27, 2002, and which claimed priority from United States Provisional 
patent application No. 60/171,606 filed December 23, 1999, Canadian patent 
application No. 2,314,723 filed July 31, 2000 and Canadian patent application No. 
10 2,319,949 filed September 20, 2000. The contents of these applications are 
incorporated herein by reference. 

BACKGROUND OF THE INVENTION 

TECHNICAL FIELD 
15 The invention relates to optical devices and is especially applicable to 

waveguide structures and integrated optics. 

BACKGROUND ART 

This specification refers to several published articles. For convenience, the 
20 articles are cited in full in a numbered list at the end of the description and cited by 
number in the specification itself. The contents of these articles are incorporated herein 
by reference and the reader is directed to them for reference. 

In the context of this patent specification, the term "optical radiation" embraces 
electromagnetic waves having wavelengths in the infrared, visible and ultraviolet 
25 ranges. 

The terms "finite" and "infinite" as used herein are used by persons skilled in 
this art to distinguish between waveguides having "finite" widths in which the actual 
width is significant to the performance of the waveguide and the physics governing its 
operation and so-called "infinite" waveguides where the width is so great that it has no 

30 significant effect upon the performance and physics or operation. 

At optical wavelengths, the electromagnetic properties of some metals closely 
resemble those of an electron gas, or equivalently of a cold plasma. Metals that 
resemble an almost ideal plasma are commonly termed "noble metals" and include, 
among others, gold, silver and copper. Numerous experiments as well as classical 

35 electron theory both yield an equivalent negative dielectric constant for many metals 
when excited by an electromagnetic wave at or near optical wavelengths [1 ,2], In a 
recent experimental study, the dielectric function of silver has been accurately 
measured over the visible optical spectrum and a very close correlation between the 
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measured dielectric function and that obtained via the electron gas model has been 
demonstrated [3]. 

It is well-known that the interface between semi-infinite materials having 
positive and negative dielectric constants can guide TM (Transverse Magnetic) surface 
5 waves. In the case of a metal-dielectric interface at optical wavelengths, these waves 
are termed plasmon-polari ton modes and propagate as electromagnetic fields coupled to 
surface plasmons (surface plasma oscillations) comprised of conduction electrons in the 
metal [4]. 

It is known to use a metal film of a certain thickness bounded by dielectrics 

10 above and below as an optical slab (planar, infinitely wide) waveguiding structure, 

with the core of the waveguide being the metal film. When the film is thin enough, the 
plasmon-polariton modes guided by the interfaces become coupled due to field 
tunnelling through the metal, thus creating supermodes that exhibit dispersion with 
metal thickness. The modes supported by infinitely wide symmetric and asymmetric 

15 metal film structures are well-known, as these structures have been studied by 

numerous researchers; some notable published works include references [4] to [10]. 

In general, only two purely bound TM modes, each having three field 
components, are guided by an infinitely wide metal film waveguide. In the plane 
perpendicular to the direction of wave propagation, the electric field of the modes is 

20 comprised of a single component, normal to the interfaces and having either a 
symmetric or asymmetric spatial distribution across the waveguide. Consequently, 
these modes are denoted s b and a b modes, respectively. The s b mode can have a small 
attenuation constant and is often termed a long-range surface plasmon-polariton. The 
fields related to the a b mode penetrate further into the metal than in the case of the s b 

25 mode and can be much lossier by comparison. Interest in the modes supported by thin 
metal films has recently intensified due to their useful application in optical 
communications devices and components. Metal films are commonly employed in 
optical polarizing devices [11] while long-range surface plasmon-polaritons can be used 
for signal transmission [7]. In addition to purely bound modes, leaky modes are also 

30 known to be supported by these structures. 

Infinitely wide metal film structures, however, are of limited practical interest 
since they offer one-dimensional (1-D) field confinement only, with confinement 
occurring along the vertical axis perpendicular to the direction of wave propagation, 
implying that modes will spread out laterally as they propagate from a point source 

35 used as the excitation. Metal films of finite width have recently been proposed in 
connection with polarizing devices [12], but merely as a cladding. 

In addition to the lack of lateral confinement, plasmon-polariton waves guided 
by a metal -dielectric interface are in general quite lossy. Even long-range surface 



3 



plasmons guided by a metal film can be lossy by comparison with dielectric 
waveguides. Known devices exploit this high loss associated with surface plasmons for 
the construction of plasmon-polariton based modulators and switches. Generally, 
known plasmon-polariton based modulator and switch devices can be classified along 
5 two distinct architectures. The first architecture is based on the phenomenon of 

attenuated total reflection (ATR) and the second architecture is based on mode coupling 
between a dielectric waveguide and a nearby metal. Both architectures depend on the 
dissipation of optical power within an interacting metal structure. 

ATR based devices depend on the coupling of an optical beam , which is 

10 incident upon a dielectric-metal structure placed in optical proximity, to a surface 
plasmon-polariton mode supported by the metal structure. At a specific angle of 
incidence, which depends on the materials used and the particular geometry of the 
device, coupling to a plasmon mode is maximised and a drop in the power reflected 
from the metal surface is observed. ATR based modulators make use of this attenuated 

15 reflection phenomenon along with means for varying electrically or otherwise at least 
one of the optical parameters of one of the dielectrics bounding the metal structure in 
order to shift the angle of incidence where maximum coupling to plasmons occurs. 
Electrically shifting the angle of maximum coupling results in a modulation of the 
intensity of the reflected light. Examples of devices that are based on this architecture 

20 are disclosed in references [23] to [36]. 

Mode coupling devices are based on the optical coupling of light propagating in 
a dielectric waveguide to a nearby metal film placed a certain distance away and in 
parallel with the dielectric waveguide. The coupling coefficient between the optical 
mode propagating in the waveguide and the plasmon-polariton mode supported by the 

25 nearby metal film is adjusted via the materials selected and the geometrical parameters 
of the device. Means is provided for varying, electrically or otherwise, at least one of 
the optical parameters of one of the dielectrics bounding the metal. Varying an optical 
parameter (the index of refraction, say) varies the coupling coefficient between the 
optical wave propagating in the dielectric waveguide and the lossy plasmon-polariton 

30 wave supported by the metal. This results in a modulation in the intensity of the light 
exiting the dielectric waveguide. References [37] to [40] disclose various device 
implementations based upon this phenomenon. Reference [41] further discusses the 
physical phenomenon underlying the operation of these devices. 

Reference [42] discusses an application of the ATR phenomenon for realising 

35 an optical switch or bistable device. 

These known modulation and switching devices disadvantageous^ require 
relative high control voltages and have limited electrical/optical bandwidth. 
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SUMMARY OF THE INVENTION 
The present invention seeks to eliminate, or at least mitigate, one or more of the 
disadvantages of the prior art, or at least provide an alternative. 

According to one aspect of the present invention there is provided an optical 
5 device comprising a waveguide structure formed by a thin strip of material having a 
relatively high free charge carrier density surrounded by material having a relatively 
low free charge carrier density, the strip having finite width and thickness with 
dimensions such that optical radiation having a free space wavelength in the range from 
about 0.8nm to about 2|im couples to the strip and propagates along the length of the 
10 strip as a plasmon-polariton wave. 

Such a strip of finite width offers two-dimensional (2-D) confinement in the 
transverse plane, i.e. perpendicular to the direction of propagation, and, since suitable 
low-loss waveguides can be fabricated from such strip, it may be useful for signal 
transmission and routing or to construct components such as couplers, power splitters, 
15 interferometers, modulators, switches and other typical components of integrated 
optics. In such devices, different sections of the strip serving different functions, in 
some cases in combination with additional electrodes. The strip sections may be 
discrete and concatenated or otherwise interrelated, or sections of one or more 
continuous strips. 

20 For example, where the optical radiation has a free-space wavelength of 1550 

nm, arid the waveguide is made of a strip of a noble metal surrounded by a good 
dielectric, say glass, suitable dimensions for the strip are thickness less than about 0. 1 
microns, preferably about 20 nm, and width of a few microns, preferably about 4 
microns. 

25 The strip could be straight, curved, bent, tapered, and so on. 

The dielectric material may be inhomogeneous, for example a combination of 
slabs, strips, laminae, and so on. The conductive or semiconductive strip may be 
inhomogeneous, for example a layer of gold and a layer of titanium. 

The plasmon-polariton wave which propagates along the structure may be 
30 excited by an appropriate optical field incident at one of the ends of the waveguide, as 
in an end-fire configuration, and/or by a different radiation coupling means. 

The device may further comprise means for varying the value of an 
electromagnetic property of at least a portion of said surrounding material so as to vary 
the propagation characteristics of the waveguide structure and the propagation of the 
35 plasmon-polariton wave. 

In some embodiments of the invention, for one said value of the 
electromagnetic property, propagation of the plasmon-polariton wave is supported and, 
for another value of said electromagnetic property, propagation of the plasmon- 
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polariton wave is at least inhibited. Such embodiments may comprise modulators or 
switches. 

Different embodiments of the invention may employ different means of varying 
the electromagnetic property, such as varying the size of at least one of said portions, 
5 especially if it comprises a fluid. 

The at least one variable electromagnetic property of the material may comprise 
permittivity , permeability or conductivity . 

Various objects, features, aspects and advantages of the present invention will 
become more apparent from the following detailed description, taken in conjunction 
10 with the accompanying drawings, of preferred embodiments of the invention which are 
described by way of example only. 

BRIEF DESCRIPTION OF DRAWINGS 

Figures 1(a) and 1(b) are a cross-sectional illustration and a plan view, 

15 respectively, of a symmetric waveguide structure embodying the present invention in 
which the core is comprised of a lossy metal film of thickness r, width u>, length / and 
permittivity e 2 embedded in a cladding or background comprising an "infinite" 
homogeneous dielectric having a permittivity e x ; 

Figures 2(a) and 2(b) illustrate dispersion characteristics with thickness of the 

20 first eight modes supported by a symmetric metal film waveguide of width w = l\im. 
The a^, and s b modes supported for the case w = « are shown for comparison, (a) 
Normalized phase constant; (b) Normalized attenuation constant; 

Figures 3(a),(b),(c),(d),(e) and (f) illustrate the spatial distribution of the six 
field components related to the ss b mode supported by a symmetric metal film 

25 waveguide of thickness / = lOOnm and width w = lfim. The waveguide cross-section 
is located in the x - y plane and the metal is bounded by the region -0.5 < x < 0.5[im 
and -0.05 < y < 0.05nm, outlined as the rectangular dashed contour. The field 
distributions are normalized such that max\Re{E y }\ = 1; 

Figures 4(a),(b),(c),(d),(e) and (f) illustrate the spatial distribution of the six 

30 field components related to the sal mode supported by a symmetric metal film 

waveguide of thickness t = \00nm and width w = l\im. The waveguide cross-section 
is located in the x - y plane and the metal is bounded by the region -0.5 < x % 0.5^ra 
and -0.05 < y < 0.05jim, outlined as the rectangular dashed contour. The field 
distributions are normalized such that mcalRe^W = 1; 

35 Figures 5(a),(b),(c),(d),(e) and (f) illustrate the spatial distribution of the six 

field components related to the a$l mode supported by a symmetric metal film 
waveguide of thickness t = lOOnm and width w = l[xm. The waveguide cross-section 
is located in the x - y plane and the metal is bounded by the region -0.5 % x < 0.5|im 
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and -0.05 < y <; O.OSura, outlined as the rectangular dashed contour. The field 
distributions are normalized such that max\Re{Ey}\ = 1; 

Figures 6(a),(b),(c),(d),(e) and (f) illustrate the spatial distribution of the six 
field components related to the aa£ mode supported by a symmetric metal film 
5 waveguide of thickness t = lOOnm and width w = \\un. The waveguide cross-section 
is located in the x - y plane and the metal is bounded by the region -0.5 % x <, 0.5pm 
and -0.05 ^ y <. 0. 05 \im, outlined as the rectangular dashed contour. The field 
distributions are normalized such that wax \Re{E y }\ = 1; 

Figures 7(a),(b),(c),(d),(e) and (f) are contour plots of Re{S z } associated with 
10 the ss£ mode for symmetric metal film waveguides of width w = l[xm and various 
thicknesses. The power confinement factor cf is also given in all cases, and is 
computed via equation (12) with the area of the waveguide core A c taken as the area of 
the metal region. In all cases, the outline of the metal film is shown as the rectangular 
dashed contour; 

15 Figure 8 illustrates a normalized profile of Re{S z } associated with 

the ss£ mode for a symmetric metal film waveguide of width w = l\im and thickness 
t = 2Qnm. The waveguide cross-section is located in the x'-y plane and the metal film 
is bounded by the region -0.5 < x < O.Sjim and -0.01 < y s O.Ol^m, outlined as the 
rectangular dashed contour; 

20 Figures 9(a),(b),(c) and (d) illustrate the spatial distribution of the field 

component related to some higher order modes supported by a symmetric metal film 
waveguide of thickness t = lOOrcm and width w = l[im. In all cases, the waveguide 
cross-section is located in the x - y plane and the metal film is bounded by the region - 
0.5 % x < 0.5|im and -0.05 <; y < 0.05\im, outlined as the rectangular dashed contour; 

25 Figures 10(a) and (b) illustrate dispersion characteristics with thickness of the 

first six modes supported by a symmetric metal film waveguide of width w = 0.5(i/n. 
The a b and s b modes supported for the case w = «> are shown for comparison, (a) 
Normalized phase constant; (b) Normalized attenuation constant; 

Figures 1 1(a) and (b) illustrate dispersion characteristics with thickness of 

30 the ss% mode supported by symmetric metal film waveguides of various widths. The 
s b mode supported for the case w = « is shown for comparison, (a) Normalized phase 
constant; (b) Normalized attenuation constant; 

Figures 12(a), (b),(c) and (d) illustrate a contour plot of Re{S z ] associated with 
the ss£ mode for symmetric metal film waveguides of thickness t = 2Qnm and various 

35 widths. The power confinement factor cf is also given in all cases, and is computed 
via equation (12) with the area of the waveguide core A c taken as the area of the metal 
region. In all cases, the outline of the metal film is shown as the rectangular dashed 
contour; 
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Figure 13 illustrates dispersion characteristics with thickness of the ss£ mode 
supported by a symmetric metal film waveguide of width w = 0.5pm for various 
background permittivities e r l . The normalized phase constant is plotted on the left axis 
and the normalized attenuation constant is plotted on the right one; 
5 Figures 14(a),(b),(c) and (d) illustrates a contour plot of Re{Sz} associated with 

the ss£ mode for a symmetric metal film waveguide of width w = 0.5pm and 
thickness t — 20nm for various background permittivities e r V In all cases, the outline 
of the metal film is shown as the rectangular dashed contour; 

Figures 15(a) and (b) illustrate dispersion characteristics with frequency of 
10 the ss£ mode supported by symmetric metal film waveguides of width w = 0.5pm 
and w = 1pm and various thicknesses t. The s b mode supported for the case w = <*> 
and the thicknesses considered is shown for comparison, (a) Normalized phase 
constant, (b) Mode power attenuation computed using Equation (16) and scaled to 
dB/cm; 

15 Figures 16(a),(b),(c),(d),(e) and (f) illustrate a contour plot of Re{S z } associated 

with the ss£ mode for symmetric metal film waveguides of width w = 0.5pm and w 
= 1pm, and thickness t — 20nm at various free-space wavelengths of excitation A 0 . In 
all cases, the outline of the metal film is shown as the rectangular dashed contour; 
Figures 17(a) and 17(b) are a cross-sectional view and a plan view, 

20 respectively, of a second embodiment of the invention in the form of an asymmetric 
waveguide structure formed by a core comprising a lossy metal film of thickness f, 
width w and permittivity e 2 supported by a homogeneous semi-infinite substrate of 
permittivity e t and with a cover or superstate comprising a homogeneous semi-infinite 
dielectric of permittivity e 3 ; 

25 Figures 18(a) and 18(b) illustrate dispersion characteristics with thickness of the 

first seven modes supported by an asymmetric metal film waveguide of width w = 1 
pm. The a b and s b modes supported for the case w = . « are shown for comparison, (a) 
Normalized phase constant, (b) Normalized attenuation constant; 

Figures 19(a), (b),(c) and (d) illustrate spatial distribution of the Ey field 

30 component related to the ss b ° mode supported by an asymmetric metal film waveguide 
of width w = 1 pm for four film thicknesses. The waveguide cross-section is located in 
the x-y plane and the metal region is outlined as the rectangular dashed contour. The 
field distributions are normalized such that max\di{E y } \ = 1; 

Figures 20(a), (b),(c) and (d) illustrate spatial distribution of the E y field 

35 component related to two higher order modes supported by an asymmetric metal film 
waveguide of width w = 1 \im for two film thicknesses. In all cases, the waveguide 
cross-section is located in the x-y plane and the metal region is outlined as the 
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rectangular dashed contour. The field distributions are normalized such that 
max|9t{Ey}| =1; 

Figures 21(a) and (b) illustrate dispersion characteristics with thickness of the 
first six modes supported by an asymmetric metal film waveguide of width w = 1 [im. 
5 The a b and s b modes supported for the case w = °° are shown for comparison, (a) 
Normalized phase constant, (b) Normalized attenuation constant; 

Figures 22(a), (b),(c) and (d) illustrate spatial distribution of the Ey field 
component related to modes supported by an asymmetric metal film waveguide of 
width w = 1 ^m. In all cases, the waveguide cross-section is located in the x-y plane 
10 and the metal region is outlined as the rectangular dashed contour. The field 
distributions are normalized such that max\^{E y } | = 1; 

Figures 23(a) and 23(b) illustrate dispersion characteristics with thickness of the 
first six modes supported by an asymmetric metal film waveguide of width w = 0.5 
[im. The a b and \ modes supported for the case w = °° are shown for comparison, (a) 
15 Normalized phase constant, (b) Normalized attenuation constant; 

Figures 24(a) and 24(b) illustrate dispersion characteristics with thickness of 
the ss£ and sal m °des supported by an asymmetric metal film waveguide of width w 
= 0.5 [im for various cases of € 3 . (a) Normalized phase constant; the inset shows an 
enlarged view of the region bounded by 0.04 < t s 0.08 [im and 2.0 <. p/p 0 <> 2.3. (b) 
20 Normalized attenuation constant; the inset shows an enlarged view of the region 
bounded by 0.05 < t s 0.08 \im and 7.0 X 10~ 3 < a/p 0 * 2.0 X 10' 2 ; 

Figures 25(a),(b),(c) and (d) illustrate spatial distribution of the Ey field 
component related to the sa b l mode supported by an asymmetric metal film waveguide 
of width w = 0.5 [im for four film thicknesses. The waveguide cross-section is located 
25 in the x-y plane and the metal region is outlined as the rectangular dashed contour. 
The field distributions are normalized such that max\ 9t{Ey} | = 1; 

Figures 26(a), (b),(c) and (d) illustrate a contour plot of St{S z } associated with 
the long-ranging modes supported by asymmetric metal film waveguides of width w = 
0.5 [im and having different superstate permittivities e 3 . In all cases, the outline of the 
30 metal film is shown as the rectangular dashed contour; 

Figure 27 is a plan view of a waveguide with opposite sides stepped to provide 
different widths; 

Figure 28 is a plan view of a waveguide which is tapered and slanted; 

Figure 29 is a plan view of a trapezoidal waveguide; 
35 Figure 30 is a plan view of a waveguide having curved side edges and suitable 

for use as a transition piece; 

Figure 31 is a plan view of a curved waveguide section suitable for 
interconnecting waveguides at a corner; 
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Figure 32 is a plan view of a two-way splitter/combiner formed by a 
combination of three straight waveguide sections and one tapered waveguide section; 

Figure 33(a) is a plan view of an angled junction using a slanted section; 

Figure 33(b) is a plan view of an offset junction using an S-bend; 
5 Figure 34(a) is a plan view of a power divider formed by a trapezoidal section 

and pairs of concatenated bends; 

Figure 34(b) is a plan view of a power divider similar to that shown in Figure 
34(a) but with a transition section formed by mirroring and overlapping curved 
sections; 

10 Figure 35 is a plan view of a Mach-Zender interferometer formed using a 

combination of the waveguide sections; 

Figure 36(a) is a schematic plan view of a modulator using the Mach-Zender 
waveguide structure of Figure 35; 

Figures 36(b) and 36(c) are inset diagrams illustrating alternative ways of 
15 applying a modulation control voltage; 

Figure 37 is a plan view of a modulator using the Mach-Zender waveguide 
structure of Figure 35 and illustrating magnetic field control; 

Figure 38 is a plan view of a periodic structure formed by a series of unit cells 
each comprising two waveguide sections having different widths and lengths; 
20 Figure 39 is a plan view of a periodic waveguide structure formed by a series of 

unit cells each comprising two opposed trapezoidal waveguide sections; 

Figure 40(a) is a plan view of an edge coupler formed by two parallel strips of 
straight waveguide with various other waveguides for coupling signals to and from 
them; 

25 Figure 40(b) is an inset diagram illustrating a way of applying a modulation 

control voltage; 

Figure 40(c) is a plan view of an edge coupler similar to that shown in Figure 
40(a) but using S-bends; 

Figure 41(a) is a perspective view of an edge coupler in which the parallel strips 
30 are not co-planar; 

Figure 41(b) is an inset diagram illustrating a way of applying a modulation 
control voltage; 

Figure 42 is a plan view of an intersection formed by four sections of 
waveguide; . 

35 Figures 43(a) and 43(b) are a schematic front view and corresponding top plan 

view of an electro-optic modulator employing the waveguide structure of Figure 17(a); 
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Figures 44(a) and 44(b) are a schematic front view and corresponding top view 
of an alternative electro-optic modulator also using the waveguide structure of Figure 
17(a); 

Figure 44(c) illustrates an alternative connection arrangement of the modulator 
5 of Figure 44(a); 

Figure 45 is a schematic front view of a third embodiment of electro-optic 
modulator also using the waveguide structure of Figure 17(a); 

Figure 46 is a schematic front view of a magneto-optic modulator also using the 
waveguide structure of Figure 17(a); 
10 Figure 47 is a schematic front view of a thermo-optic modulator also using the 

waveguide structure of Figure 17(a); 

Figure 48 is a schematic perspective view of an electro-optic switch also using 
the waveguide structure of Figure 17(a); 

Figure 49 is a schematic perspective view of a magneto-optic switch also using 
15 the waveguide structure of Figure 17(a); 

Figure 50 is a schematic perspective view of a thermo-optic switch also using 
the waveguide structure of Figure 17(a); 

Figure 51 gives the mode power attenuation for metal film waveguides of 
various widths and thicknesses. The metal used is Au and the background dielectric is 
20 Si0 2 . The optical free-space wavelength of analysis is set to X 0 = 1.55 \xm; and 

Figure 52 gives the mode power attenuation for metal Film waveguides of 
various widths and thicknesses. The metal used is Al and the background dielectric is 
Si0 2 . The optical free-space wavelength of analysis is set to X 0 = 1.55 nm. 

25 DESCRIPTION OF THE PREFERRED EMBODIMENT 
I. Introduction 

In order to facilitate an understanding of the specific optical devices embodying 
the invention, their theoretical basis will first be explained with reference to Figures 1 
to 26(d). 

30 The following is a comprehensive description of the purely bound modes of 

propagation supported by symmetric and asymmetric waveguide structures comprised 
of a thin lossy metal film of /m/Ye-width as the core. The core can be embedded in an 
"infinite" homogeneous dielectric medium as shown in Figure 1(a) or supported by a 
semi-infinite homogeneous dielectric substrate and covered by a different semi-infinite 

35 homogeneous dielectric superstate as shown in Figure 17(a). The description is 
organized as follows. Section II summarizes the physical basis and numerical 
technique used to analyze the structures of interest. Sections III through VI describe 
the modes supported by symmetric structures as shown in Figure 1(a) and sections VII 
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through X describe the modes supported by asymmetric structures as shown in Figure 
17(a). Concluding remarks are given in section XI. 

II. Physical. Basis and Numerical Technique 
5 A symmetric structure embodying the present invention is shown in Figures 1(a) 

and 1(b). It comprises a lossy metal film of thickness t, width w and equivalent 
permittivity e 2 , surrounded by a cladding or background comprising an infinite 
homogeneous dielectric of permittivity e,. Figure 17(a) shows an asymmetric structure 
(z x * e 3 ) embodying the present invention. The Cartesian coordinate axes used for the 
10 analysis are also shown with propagation taking place along the z axis, which is out of 
the page. 

It is assumed that the metal region shown in Figures 1(a) and 17(a) can be 
modelled as an electron gas over the wavelengths of interest. According to classical or 
Drude electron theory, the complex relative permittivity of the metal region is given by 
15 the well-known plasma frequency dispersion relation [4]: 



8 r,2 



( 2 \ 

i w ' 



0> 2 + V 2 



i 2 ^ 

k o>(a> 2 + v 2 



(1) 



where a> is the excitation frequency, co p is the electron plasma frequency and v is the 
20 effective electron collision frequency, often expressed as v = 1/x with x defined as 
the relaxation time of electrons in the metal. When o> 2 +v 2 < to* (which is the case 
for many metals at optical wavelengths) a negative value for the real part e r 2 is 
obtained, implying that plasmon-polariton modes can be supported at interfaces with 
normal dielectrics. 

25 ' 

Electromagnetic Wave and Field Equations 

The modes supported by the structures are obtained by solving a suitably 
defined boundary value problem based on Maxwell's equations written in the frequency 
domain for a lossy inhomogeneous isotropic medium! Uncoupling Maxwell's 
30 equations yields the following time-harmonic vectorial wave equations for the E and H 
fields: 

v x v x E - o) 2 e(x,y)(iE = 0 (2) 
v x e(x,y)- 1 7 x H - g>VH = 0 (3) 
where the permittivity e is a complex function of cross-sectional space, and describes 
35 the waveguide structure. For the structures analyzed in this description, ^ is 
homogeneous and taken as the permeability of free space |i 0 . 

Due to the nature of the numerical method used to solve the boundary value 
problem, the implicit y dependence of the permittivity can be immediately removed 
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since any inhomogeneity along y is treated by dividing the structure into a number of 
layers that are homogeneous along this direction, and suitable boundary conditions are 
applied between them. 

The two vectorial wave equations (2) and (3) are expanded in each layer: into 
5 scalar wave equations, some being coupled by virtue of the remaining inhomogeneity 
in e along x. Since the structure under consideration is invariant along the propagation 
axis (taken to be in the +z direction), the mode fields vary along this dimension 
according to e ~ YZ where y = a + is the complex propagation constant of the 
mode, a being its attenuation constant and p its phase constant. Substituting this field 
10 dependency into the scalar wave equations, and writing them for TE X (E x = 0) and 
TM X (H x = 0) modes while making use of v*[e(x)E] = 0 and v H = 0 accordingly, 
yields simplified and uncoupled scalar wave equations that are readily solved. The Ey 
component of the TE X modes must satisfy the Helmholtz wave equation: 
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d -E™ + JL E? + [Y 2 + « V»] Ey™ = 0 



dx 



dy 



(4) 
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and the H y component of the TM X modes must satisfy the Sturm-Liouville wave 
equation: 
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dx 


e(x) dx 
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d -H™ + tf = «?»e{x)\H y m = Q 



dy 



(5) 



The superposition of the TE X and TM X mode families then describes any mode 
propagating in the structure analyzed. The electric and magnetic field components 

resulting from this superposition are given by the following equations: 

a tt tm) 2 tt ™ 
— H y I +bi z \iH y 



E = 



-1 
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dx\ e(x) dx 



(6) 
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(7) 
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E - ±±E™ + -J— ±H y m 
y dy ywe(x) dx 



(8) 
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H = 
y j<*>Y\i dxdy 



E* * H™ 



J(*>\l dx 



1 d c te 1 d rj m 
E + —-H 

Y dy 



(9) 
(10) 
(11) 
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In order to obtain a mode of propagation supported by a waveguiding structure, 
the Helmholtz and Sturm-Liouville wave equations (4) and (5), along with the field 
equations (6)-(ll), must be solved for the propagation constant y using appropriate 
boundary conditions applied between layers and at the horizontal and vertical limits. 

5 

Poynting Vector and Power Confinement Factor 

The power confinement factor is defined as the ratio of mode complex power 
carried through a portion of a waveguide's cross-section with respect to the mode 
complex power carried through the entire waveguide cross-section. . Formally it is 
10 expressed as: 



20 



cf = 



Sds\ 



(12) 



Ik* 

where A c is usually taken as the area of the waveguide core and A„ implies integration 
15 over the entire waveguide cross-section (which can be all cross- sectional space for an 
open structure) or the entire cross-sectional computational domain. S z refers to the z 
component of the Poynting vector: 



and H* y denotes the complex conjugate of H x y . The spatial distribution of a 
component of the Poynting vector is easily computed from the spatial distribution of 
the relevant electric and magnetic mode field components. 

25 Numerical Solution Approach 

The boundary value problem governed by equations (4) to (11) is solved by 
applying the Method of Lines (MoL). The MoL is a well-known numerical technique 
and its application to various electromagnetic problems, including optical waveguiding, 
is well-established [14]. The MoL is rigorous, accurate and flexible. It can handle a 

30 wide variety of waveguide geometries, including the structures at hand. The method is 
not known to generate spurious or non-physical modes. The MoL formulation used 
herein is based on the formulation reported in [15], but simplified for isotropic media, 
as prescribed by equations (4) - (11) and reported in [16]. Except for a 1-D spatial 
discretization, the method is exact. 

35 The main idea behind the MoL is that the differential field equations governing 

a waveguiding problem are discretized only as far as necessary so that generalized 
analytic solutions can be applied to derive a homogeneous matrix problem describing 
all modes supported by the structure. This approach renders the method accurate and 
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computationally efficient since only -1 dimensions must be discretized to solve an N 
dimension problem. In the case of a two-dimensional (2-D) waveguiding structure, 
this means that only one spatial dimension needs to be discretized. The main features 
of this procedure, as applied to a modal analysis problem, are described below. 
5 • The x axis and the function e(x) are discretized using two shifted non- 

equidistant line systems, parallel to the y axis. 

• The differential operators d\dx and d 2 /d 2 x in the wave and field equations 
are replaced by finite difference approximations that include the lateral 
boundary conditions. 

10 • The discretized wave equations are diagonalized using appropriate 

transformations matrices. 

• The diagonalization procedure yields in the transform domain two 
systems of uncoupled one-dimensional (1-D) differential equations along 
the remaining dimension (in this case along the y axis). 

15 • These differential equations are solved analytically and tangential field 

matching conditions are applied at interfaces between layers along with 
the top and bottom boundary conditions. 

• The last field matching condition, applied near the center of the 
structure, yields a homogeneous matrix equation of the 

20 , form G(y)i = 0 which operates on transformed tangential fields. 

• The complex propagation constant y of modes is then obtained by 
searching for values that satisfy det[G(y)] = 0. 

• Once the propagation constant of a mode has been determined, the 
spatial distribution of all six field components of the mode are easily 

25 generated. 

• A mode power confinement factor can be computed by first computing 
the spatial distribution of S z which is then integrated according to 
Equation (12). 

The open structures shown in Figures 1(a) and 17(a) are discretized along the* 
30 axis and the generalized analytic solution applied along the y axis. The physical 
symmetry of the structures is exploited to increase the accuracy of the results and to 
reduce the numerical effort required to generate the mode solutions. For the 
symmetric structure shown in Figure 1(a), this is achieved by placing either electric 
wall (E^ = 0) or magnetic wall (H^ = 0) boundary conditions along the x and y axes. 
35 For the asymmetric structure shown in Figure 17(a), this is achieved by placing electric 
wall or magnetic wall boundary conditions along the y axis only. The remaining 
horizontal boundary conditions are placed at infinity and the remaining lateral 
boundary condition is either placed far enough from the guide to have a negligible 
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effect on the mode calculation, or a lateral absorbing boundary condition is used to 
simulate infinite space, depending on the level of confinement observed in the resulting 
mode. The use of numerical methods to solve differential equations inevitably raises 
questions regarding the convergence of computed results and their accuracy. The 
5 propagation constant of a mode computed using the method of lines converges in a 
monotonic or smooth manner with a reduction in the discretization interval (which 
increases the number of lines in the calculation and thus the numerical effort). This 
suggests that extrapolation can be used to generate a more accurate value for the 
propagation constant, and this value can then be used to compute the error in values 

10 obtained using the coarser discretizations [17]. This anticipated error does not 
correspond to the actual error in the propagation constant as the latter could only be 
known if the analytic or exact value were available. The anticipated error, however, 
still provides a useful measure of accuracy since it must tend toward zero as more 
accurate results are generated. 

15 The convergence of the computed propagation constant of the modes supported 

by the structures of interest has been monitored during the entire study. The 
anticipated error in the results presented herein is estimated as 1 % on average and 6% 
in the worst case. These error values are based on extrapolated propagation constants 
computed using Richardson's extrapolation formula [18]. 

20 

III. Mode Characteristics and Evolution With Film Thickness: Symmetric Structures 
A. Review of Mode Solutions for Metal Film Slab Waveguides 

The review begins with the reproduction of results for an infinitely wide 
symmetric metal film waveguide, as shown in Figure 1(a) with w = ~ 5 taken from the 
25 standard work on such structures [6]. In order to remain consistent with their results, 
the optical free-space wavelength of excitation is set to X 0 — 0.633^m and their value 
for the relative permittivity of the silver film at this wavelength is used: e r 2 = -19 - 
j0.53. The relative permittivity of the top and bottom dielectric regions is set to e r l = 
4. 

30 An infinitely wide structure supports only two purely bound TM (E x = H y = 

H z = 0) modes having transverse field components Ey and H x that exhibit asymmetry 
or symmetry with respect to the x axis. These modes are created from the coupling of 
individual plasmon-polariton modes supported by the top and bottom interfaces and 
they exhibit dispersion with film thickness. The widely accepted nomenclature for 

35 identifying them consists in using the letters a or s for asymmetric or symmetric 

transverse field distributions, respectively, followed by a subscript b or / for bound or 
leaky modes, respectively. The propagation constants of the a b and s b modes have been 
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computed as a function of film thickness and the normalized phase and attenuation 
constants are plotted in Figures 2(a) and 2(b), respectively. 

From Figures 2(a) and 2(b), it is observed that the a b and s b modes become 
degenerate with increasing film thickness. As the separation between the top and 
bottom interfaces increases, the a b and s b modes begin to split into a pair of uncoupled 
plasmon-polariton modes localized at the metal-dielectric interfaces. The propagation 
constants of the a b and s b modes thus tend towards that of a plasmon-polariton mode 
supported by the interface between semi-infinite metallic and dielectric regions, which 
is given via the following equations [6]: 



10 



P/P 0 = -X*\ 
a/p 0 = -Im 



6 M € r,2 



€ r, 1 + 6 r ,7 



6 M € r,2 



r,l r,2 



(14) 
(15) 



15 where p 0 = co/c 0 with c 0 being the velocity of light in free space, and e r l and e r 2 are 
the complex relative permittivities of the materials used. Using the above equations, 
values of p/p 0 = 2.250646 and <x/p o = 0.836247 x 10 2 are obtained for e r l = 4 and 
e r2 = -19-j0.53. 

As the thickness of the film decreases, the phase and attenuation constants of 
20 the a b mode increase, becoming very large for very thin films. This is due to the fact 
that the fields of this mode penetrate progressively deeper into the metal as its thickness 
is reduced. In the case of the s b mode, a decreasing film thickness causes the opposite 
effect, that is, the fields penetrate progressively more into the top and bottom dielectric 
regions and less into the metal. The propagation constant of this mode thus tends 
25 asymptotically towards that of a TEM (Transverse ElectroMagnetic) wave propagating 
in an infinite medium having the same permittivity as the top and bottom dielectric 
regions. In this case, the attenuation constant decreases asymptotically towards zero 
since losses were neglected in these regions. The a b and s b modes do not have a cutoff 
thickness. 

30 The fields in an infinitely wide structure do not exhibit any spatial variation 

along jc. Due to the nature of the MoL, and to the fact that the generalized analytical 
solution is applied along they dimension, our results do not contain discretization 
errors and thus are in perfect agreement with those reported in [6], 

35 B. Modes Supported by a Metal Film of Width w = l|im 

Next, the analysis of the structure shown in Figure 1(a) for the case w = 1pm 
will be explained. The material parameters and free-space wavelength that were used 
in the previous case (w — <*>) were also used here. The MoL was applied and the 
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discretization adjusted until convergence of the propagation constant was observed. 
The physical quarter-symmetry of the structure was exploited by placing vertical and 
horizontal electric or magnetic walls along they and x axes, respectively, which leads 
to four possible wall combinations as listed in Table 1 . The first two purely bound 
5 (non-leaky) modes for each wall combination were found and their dispersion with 
metal thickness computed. The results for these eight modes are shown in Figures 2(a) 
and 2(b). 

Table 1: Vertical-Horizontal wall combinations used along the axes of symmetry and 
10 proposed mode nomenclature: ew - electric wall, mw - magnetic wall. 



V-H Walls 


Mode 


ew-ew 


m 

as„ 


mw-ew 




mw-mw 


< 


ew-mw 


< 



Unlike its slab counterpart, pure TM modes are not supported by a metal film 
of finite width: all six field components are present in all modes. For a symmetric 
structure having an aspect ratio w/t > 1, the Ey field component dominates. The E x 

20 field component increases in magnitude with increasing film thickness and if w/t < 1, 
then E x dominates. It is proposed to identify the modes supported by a metal film of 
finite width, by extending the nomenclature used for metal film slab waveguides. First 
a pair of letters being a or s identify whether the main transverse electric field 
component is asymmetric or symmetric with respect to the y and x axes, respectively 

25 (in most practical structures w/t > > 1 and Ey is the main transverse electric field 

component). A superscript is then used to track the number of extrema observed in the 
spatial distribution of this field component along the largest dimension (usually along 
the x axis) between the corners. A second superscript n could be added to track the 
extrema along the other dimension (the y axis) if modes exhibiting them are found. 

30 Finally, a subscript b or / is used to identify whether the mode is bound or leaky. 
Leaky modes are known to exist in metal film slab structures and though a search for 
them has not been made at this time, their existence is anticipated. Table 1 relates the 
proposed mode nomenclature to the corresponding vertical and horizontal wall 
combinations used along the axes of symmetry. 

35 The ss£> sa^ as£ and aal modes are the first modes generated (one for each 

of the four possible quarter-symmetries listed in Table 1 , and having the largest phase 
constant) and thus may be considered as the fundamental modes supported by the 
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structure. Figures 3 to 6 show the field distributions of these modes over the cross- 
section of the waveguide for a metal film of thickness / = lOOwn. As is observed 
from these figures, the main transverse electric. field component is the E y component 
and the symmetries in the spatial distribution of this component are reflected in the 
5 mode nomenclature. The outline of the metal is clearly seen in the distribution of the 
Ey component on all of these plots. As is observed from the figures, very little field 
tunnels through the metal to couple parallel edges for this case of film thickness and 
width (very little coupling through the metal between the top and bottom edges and 
between the left and right edges), though coupling does occur along all edges between 

10 adjacent corners (mostly along the left and right ones), and also between perpendicular 
edges through the corner. 

Figures 2(a) and 2(b) suggests that the dispersion curves for these first four 
modes converge with increasing film thickness toward the propagation constant of a 
plasmon-polariton mode supported by an isolated corner (though pairs of corners in this 

15 case remain weakly coupled along the top and bottom edges due to the finite width of 
the film, even if its thickness goes to infinity). If both the film thickness and width 
were to increase further, the four fundamental modes would approach degeneracy with 
their propagation constant tending towards that of a plasmon-polariton mode supported 
by an isolated corner, and their mode fields becoming more localized near the corners 

20 of the structure with maxima occurring at all four corners and fields decaying in an 
exponential-like manner in all directions away from the corners. This is further 
supported by considering the evolution of the field distributions given in Figures 3 to 6 
as both the thickness and width increase. 

As the thickness of the film decreases, coupling between the top and bottom 

25 edges increases and the four modes split into a pair as the upper branch 

(modes sa£ and aa£ which have a dominant Ey field component exhibiting 
asymmetry with respect to the x axis) and a pair as the lower branch 
(modes ssl and as% which have a dominant Ey field component exhibiting symmetry 
with respect to the x axis), as shown in Figures 2(a) and 2(b). The pair on the upper 

30 branch remain approximately degenerate for all film thicknesses, though decreasing the 
film width would eventually break this degeneracy . The upper branch modes do not 
change in character as the film thickness decreases. Their field distributions remain 
essentially unchanged from those shown in Figures 4 and 6 with the exception that 
confinement to the metal region is increased thus causing an increase in their 

35 attenuation constant. This field behaviour is consistent with that of the a b mode 
supported by a metal film slab waveguide. 

The modes on the lower branch begin to split at a film thickness of about 80/wn, 
as shown in Figures 2(a) and 2(b). As the film thickness decreases further 
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the ss£ mode follows closely the phase and attenuation curves of the s b mode 
supported by the metal film slab waveguide. In addition to exhibiting dispersion, the 
lower branch modes change in character with decreasing thickness, their fields 
evolving from being concentrated near the corners, to having Gaussian-like 
5 distributions along the waveguide width. The E y field component of the ss£ mode 
develops an extremum near the center of the top and bottom interfaces, while that of 
the as% mode develops two extrema, one on either side of the center. Since these 
modes change in character, they should be identified when the film is fairly thick. 
Figures 7(a) to 7(f) show the evolution of the ss£ mode fields with film 

10 thickness via contour plots of /?e{S z }. S z is computed from the ss£ mode fields using 
Equation 13 and corresponds to the complex power density carried by the mode. The 
power confinement factor cf is also given in the figure for all cases, and is computed 
via equation (12) with the area of the waveguide core A c taken as the area of the metal 
region. Figures 7(a) to 7(f) clearly show how the mode fields evolve from being 

15 confined to the corners of thick films to being distributed in a Gaussian-like manner 
laterally along the top and bottom edges, as the field coupling between these edges 
increases due to a reduction in film thickness. The confinement factor becomes smaller 
as the film thickness decreases, ranging from 14% confinement to 1.6% as the 
thickness goes from 80wn to 20wn. This implies that fields become less confined to 

20 the metal, spreading out not only along the vertical dimension but along the horizontal 
one as well, as is observed by comparing Figures 7(a) and 7(b). This reduction in 
confinement to the lossy metal region explains the reduction in the attenuation constant 
of the mode with decreasing film thickness, as shown in Figure 2(b). An examination 
of all field components related to the ss£ mode reveals that the magnitudes of the 

25 weak transverse (E x , H y ) and longitudinal (E z , components decrease with 

decreasing film thickness, implying that the mode is evolving towards a TEM mode 
comprised of the E y and H x field components. Indeed, the normalized propagation 
constant of the ss% mode tends asymptotically towards the value of the normalized 
propagation constant of a TEM wave propagating in the background material (e r x = 4 

30 with no losses in this case), further supporting this fact. This field behaviour is also 
consistent with that of the s b mode supported by a metal film slab waveguide. 

Figure 8 shows the profile of Re{S z } of the ss% mode over the cross-section of 
the guide for the case t = 20/?m, providing a different perspective of the same 
information plotted as contours in Figure 7(f). Figure 8 shows that Re{S z } is negative 

35 in the metal film, implying that the mode real power is flowing in the direction 

opposite to the direction of mode propagation (or to the direction of phase velocity) in 
this region. It is clear however that the overall or net mode real power is flowing 
along the direction of propagation. It is possible that the net mode real power can be 
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made to flow in the direction opposite to that of phase velocity (as in metal film slab 
waveguides [10]) for values of e r l in the neighbourhood or greater than \Re{e r2 }\. 

Unlike the metal film slab waveguide, a metal film of finite width can support a 
number of higher order modes. The dispersion curves of the first four higher order 
5 modes (each generated from one of the symmetries listed in Table 1 are shown in 
Figures 2(a) and 2(b), and the spatial distribution of their main transverse electric field 
component is shown in Figure 9 for a film of thickness t = lOOwn. As is observed 
from Figures 9(a) to 9(d), the symmetries and number of extrema in the distributions 
of Re{E y } are reflected in the mode nomenclature. It should be noted that the nature of 

10 the nomenclature is such that all higher order modes sa b and ss b have an odd m 
while all higher order modes aa b and as b have an even ra. Comparing Figures 
9(a) to 9(d) with Figures 3(c), 4(c), 5(c) and 6(c), respectively, (ie: comparing the Ey 
component of the ssl mode shown in Figure 9(a). with the E y component of 
the ss% m °de shown in Figure 3(c), etc..) reveals that the fields of a higher order 

15 mode are comprised of the fields of the corresponding m = 0 mode with additional 
spatial oscillations or variations along the top and bottom edges of the structure due to 
the latter 1 s limited width. Making this comparison for all of the field components of 
the higher order modes found reveals this fact to be true, except for the H y field 
component which remains in all cases essentially identical to that of the corresponding 

20 m ~ 0 mode; ie: the H y field component never exhibits oscillations along the width of 
the structure. 

The evolution of the sal anc * modes with film thickness is similar to the 
evolution of the sal an( * &&1 m °des (and the a b mode supported by the metal film 
slab waveguide), in that their mode fields become more tightly confined to the metal as 
25 the thickness of the latter decreases, thereby causing an increase in the attenuation of 

1 2 

the modes, as shown in Figure 2(b). Furthermore, the sa b and aa b modes do not 
change in character with film thickness, their field distributions remaining essentially 
unchanged in appearance from those computed at a thickness of lOOwn. 

The ss b and as b modes evolve with thickness in a manner similar to the 
30 corresponding m — 0 modes (and the s b mode of the metal film slab waveguide) in the 
sense that their fields become less confined to the metal region as the thickness of the 
latter decreases, thereby reducing the attenuation of the modes as shown in Figure 2(b). 
As the thickness of the film decreases, the ss b and asl modes change in character in 
a manner similar to the corresponding m = 0 modes, their field components evolving 
35 extra variations along the top and bottom edges. 

As the thickness of the film increases, the propagation constants of the 
sal anc * ss l modes converge to a single complex value as shown in Figures 2(a) and 
2(b). This is the propagation constant of uncoupled higher order modes supported by 
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the top and bottom edges of the film. A similar observation holds for 
the aal and asl modes. The nature of these v edge modes' is clear by considering 
the evolution with increasing film thickness of the distributions shown in Figures 9(a) 
to 9(d), As the thickness of the film tends to infinity, the top edge becomes uncoupled 
5 from the bottom edge, forcing the ssl mode to become degenerate with the 

sal m °de since both have an E y field component that is symmetric with respect to the 
y axis and one extremum in its distribution along the top or bottom edge. A similar 
reasoning explains why the asl ™°de must become degenerate with the ssl m °de. 
In general, it is expected that the higher order sa™ and ss™ mode families will form 

10 degenerate pairs for a given m, as will the higher order as™ and aa™ mode 
families, with increasing film thickness. 

The aa™ and sa™ mode families do not have mode cutoff thicknesses. This 
is due to the fact that their confinement to the metal film increases with decreasing film 
thickness; thus the modes remain guided as /- 0. The as™ and ss™ mode families 

15 have cutoff thicknesses for all modes except the ssl mode, which remains guided as t 
- 0, since it evolves into the TEM mode supported by the background. The other 
modes of these families, including the asl m °de cannot propagate as t - 0 because 
their mode fields do not evolve into a TEM mode. Rather, the modes maintain 
extrema in their field distributions and such variations cannot be enforced by an infinite 

20 homogeneous medium. 

In general, the purely bound modes supported by a metal film of finite width 
appear to be formed from a coupling of modes supported by each metal-dielectric 
interface defining the structure. In a metal film of finite width, straight interfaces of 
finite length (top, bottom, left and right edges) and corner interfaces are present. Since 

25 a straight metal-dielectric interface of infinite length can support a bound plasmon- 
polariton mode then so should an isolated corner interface and a straight interface of 
finite length bounded by corners (say the edge defined by a metal of finite width 
having an infinite thickness). A preliminary analysis of an isolated corner has revealed 
that a plasmon-polariton mode is indeed supported and that the phase and attenuation 

30 constants of this mode are greater than those of the mode guided by the corresponding 
infinite straight interface, as given by Equations (14) and (15). This is due to the fact 
that fields penetrate more deeply into the metal near the corner, to couple neighbouring 
perpendicular edges. All six field components are present in such a mode, having their 
maximum value at the corner and decreasing in an exponential-like manner in all 

35 directions away from the corner. A straight interface of finite length bounded by 
corners should support a discrete spectrum of plasmon-polariton modes with the 
defining feature in the mode fields being the number of extrema in their spatial 
distribution along the edge. A mode supported by a metal film of finite width; may 
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therefore be seen as being comprised of coupled N corner modes' and s finite length 
edge modes 1 . 

The ss£ mode could be used for optical signal transmission over short 
distances. Its losses decrease with decreasing film thickness in a manner similar to the 
5 s 5 mode supported by the metal film slab waveguide. In a symmetric waveguide 
structure such as the one studied here, the ss£ mode does not have a cut-off thickness 
so losses could be made small enough to render it long-ranging, though a trade-off 
against confinement is necessary. In addition, when the metal is thin, the Ey field 
component of the mode has a maximum near the center of the metal-dielectric 
10 interfaces, with a symmetric profile similar to that shown in Figure 8. This suggests 
that the mode should be excitable using a simple end-fire technique similar to the one 
employed to excite surface plasmon-polariton modes [19,6]; this technique is based on 
maximizing the overlap between the incident field and that of the mode to be excited. 
In reference [22], the present inventor et aL disclosed that plasmon-polariton 
15 waves supported by thin metal films of finite width have recently been observed 

experimentally at optical communications wavelengths using this method of excitation. 

IV. Mode Dispersion With Film Width: Symmetric Structures 

Since the modes supported by a metal film waveguide exhibit dispersion with 
20 film thickness, it is expected that they also exhibit dispersion with film width. 

A. Modes Supported by a Metal Film of Width w = 0.5|ira 

The analysis of a metal film waveguide of width w = 0.5|im will now be 
discussed, using the material parameters and free- space wavelength that were used in 

25 the previous section. A film width of 0.5(iw was selected in order to determine the 
impact of a narrowing film on the modes supported and to demonstrate that the 
structure can still function as a waveguide though the free-space optical wavelength is 
greater than both the width and thickness of the film. 

As in the previous section, the first eight modes supported by the structure (two 

30 for each symmetry listed in Table 1) were sought, but in this case only six modes were 
found. The dispersion curves with thickness of the modes found are plotted in Figures 
10(a) and 10(b). The observations made in the previous section regarding the general 
behaviour of the modes hold for other film widths, including this one. 

2 2 

The aa b and as b modes, which were the highest order modes found for a 
35 film of width w = lpm, were not found in this case suggesting that the higher order 
modes (m > 0) in general have a cut-off width. Comparing Figure 10(a) with Figure 
2(a), it is apparent that decreasing the film width causes a decrease in the phase 
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constant of the ssl and sal socles, further supporting the existence of a cut-off 
width for these modes. 

Comparing Figures 10(a) and 10(b) with Figures 2(a) and 2(b), it is noted that 
the modes which do exhibit cutoff thicknesses (the ss™ modes with m > 0 and the 
5 as™ modes with m z 0), exhibit them at a larger thickness for a narrower film width. 
This makes it possible to design a waveguide supporting only one long-ranging mode 
(the ss£ mode) by carefully selecting the film width and thickness. 

B. Dispersion of the ss£ Mode With Film Width 

10 The dispersion with thickness of the ss£ mode is shown in Figures 11(a) and 

11(b) for numerous film widths in the range 0.25 < w < ljim, illustrating the amount of 
dispersion in the mode properties that can be expected due to a varying film width. In 
all cases, the ss£ mode evolves with decreasing film thickness into the TEM wave 
supported by the background, but this evolution occurs more rapidly for a narrower 

15 width. For a film of thickness t = 20//ra, for example, from Figure 11(a), the 

normalized phase constant of the mode supported by a film of width w = \\im is about 
2.05, while that of the mode supported by a film of width w = 0.25|im is already 
about 2. This fact is also supported by the results plotted in Figure 1 1(b) since the 
attenuation constant of the mode at a thickness of t = 20wn is closer to zero (the 

20 attenuation constant of the background) for narrow film widths compared to wider 
ones. Indeed, at a thickness of 10/zm, the attenuation of the mode for a width of w 
=0.25|im is more than an order of magnitude less than its attenuation at a width of w 
" = l\im (and more than an order of magnitude less than that of the s b mode supported 
by a metal film slab waveguide), indicating that this mode can be made even more 

25 long-ranging by reducing both the film thickness and its width. 

The dispersion of the mode with increasing film thickness also changes as a 
function of film width, as seen from Figure 11(a). This is due to the fact that the 
amount of coupling between corners along the top and bottom edges increases as the 
film narrows, implying that the mode does not evolve with increasing thickness 

30 towards a plasmon-polariton mode supported by an isolated corner, but rather towards 
a plasmon-polariton mode supported by the pair of corners coupled via these edges. 

Figures 12(a) to 12(d) show contour plots of /?<?{S Z } related to the ss£ mode 
supported by films of thickness t = 20nm and various widths. The power confinement 
factor is also given for all cases, with the area of the waveguide core A c taken as the 

35 area of the metal region. Figures 12(a) to 12(d) clearly illustrate how the fields 

become less confined to the lossy metal as its width decreases, explaining the reduction 
in attenuation shown in Figure 11(b) at this thickness. In addition, the confinement 
factor ranges from 1.64% to 0.707% for the widths considered, further corroborating 
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this fact. The fields are also seen to spread out farther, not only along the horizontal 
dimension but along the vertical one as well, as the film narrows. This indicates that 
the mode supported by a narrow film is farther along in its evolution into the TEM 
mode supported by the background, compared to a wider film of the same thickness. It 
5 is also clear from Figures 12(a) to 12(d) that the trade-off between mode confinement 
and attenuation must be made by considering not only the film thickness but its width 
as well. 

V. Effects Caused by Varying the Background Permittivity: Symmetric Structures 
10 In this section, the changes in the propagation characteristics of the ss% mode 

due to variations in the background permittivity of the waveguide. Only 
the ss£ mode is considered since the main effects are in general applicable to all 
modes. In order to isolate the effects caused by varying the background permittivity, 
the width of the metal film was fixed to w = 0.5\xm and its permittivity as well as the 
15 optical free-space wavelength of analysis were set to the values used in the previous 
sections. The relative permittivity of the background e r , was taken as the variable 
parameter. 

The dispersion with thickness of the ss£ mode is shown in Figure 13 for some 
background permittivities in the range 1 s e r , < 4. Figures 14(a) to 14(d) compare 

20 contour plots of Re{S z } related to this mode for a film of thickness t = 20/zm and for 
the same set of background permittivities used to generate the curves plotted in Figure 
13. From Figures 14(a) to 14(d); it is observed that reducing the value of the 
background permittivity causes a reduction in field confinement to the metal. This 
reduction in field confinement within the lossy metal in turn causes a reduction in the 

25 attenuation of the mode that can be quite significant, Figure 13 showing a reduction of 
almost four orders of magnitude at a film thickness of t = 20wn, as the background 
relative permittivity ranges from e r x = 4 to 1. It is also noted that the mode exhibits 
less dispersion with thickness as the background relative permittivity is reduced, since 
the normalized phase constant curves shown in Figure 13 flatten out with a reduction in 

30 the value of this parameter. 

From Figures 14(a) to 14(d), it is seen that the mode power is confined to 
within approximately one free-space wavelength in all directions away from the film in 
all cases except that shown in Figure 14(d), where fields are significant up to about two 
free-space wavelengths. In Figure 14(c), the background permittivity is roughly that of 

35 glass and from Figure 13 the corresponding normalized attenuation constant of the 
mode is about a/p 0 = 6.0 x 10" 5 . The associated mode power attenuation in dB/mm, 
computed using the following formula: 



(16) 



is about 5 dB/mm. This value of attenuation is low enough and field confinement is 
high enough as shown in Figure 14(c), to render this particular structure practical at 
5 this free-space wavelength for applications requiring short propagation lengths. 

The changes in mode properties caused by varying the background permittivity 
as discussed above are consistent with the changes observed for the modes supported 
by a metal film slab waveguide and the observations are in general applicable to the 
other modes supported by a metal film of finite width. In the case of the higher order 
10 modes (m > 0) and those exhibiting a cutoff thickness (the as™ modes for all m and 
the ss ™ modes for m > 0) additional changes in the mode properties occur. In 
particular, as the background permittivity is reduced, the cut-off widths of the higher 
order modes increase as do all relevant cut-off thicknesses. 

15 VI. Frequency Dependency of the ss% Mode Solutions: Symmetric Structures 

In order to isolate the frequency dependency of the ss£ m °de solutions, the 
geometry of the metal film was held constant and the background relative permittivity 
was set to e T l = 4. The relative permittivity of the metal film e r 2 was assumed to vary 
with the frequency of excitation according to Equation (1). In order to remain 

20 consistent with [6], the values w p = 1.29 x 10 16 rad/s and 1/v = x == 1.25 x 10 _1 V 
were adopted, though the latter do not generate exactly e r 2 = -19 - j0.53 at X 0 = 
0.633(im, which is the value used in the previous sections. This is due to the fact that 
values of o> p and t are often deduced by fitting Equation (1) to measurements. The 
values used, however, are in good agreement with recent measurements made for silver 

25 [3] and are expected to generate frequency dependent results that are realistic and 
experimentally verifiable. ' 

The dispersion characteristics of the ss% mode supported by films of width w 
= 0.5|iw and w = ljim, and thicknesses in the range 10 < t < 50nm are shown in 
Figures 15(a) and 15(b) for frequencies covering the free-space wavelength range 0.5 < 

30 X 0 z 2 \im. Curves for the s b mode supported by metal film slab waveguides (w = °°) 
of the same thicknesses are also shown for comparison. 

The results given in Figure 15(a) show that, in all cases, the normalized phase 
constant of the modes tends asymptotically towards that of the TEM wave supported by 
the background as the wavelength increases, and that the convergence to this value is 

35 steeper as the width of the film decreases (for a given thickness). The curves remain 
essentially unchanged in character as the thickness changes, but they shift upwards 
toward the top left of the graph with increasing thickness, as shown. Convergence to 
the asymptote value with increasing wavelength suggests that the ss b mode evolves 
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into the TEM mode supported by the background. It is noteworthy that the ss% mode 
can exhibit very little dispersion over a wide bandwidth, depending on the thickness 
and width of the film, though flat dispersion is also associated with low field 
confinement to the metal film. 

5 The results plotted in Figure 15(b) show in all cases a decreasing attenuation 

with increasing wavelength and the curves show a sharper drop for a narrow film (w = 
0.5\im) compared to a wide one (w = ~). The attenuation curves look essentially the 
same for all of the film thicknesses considered, though the range of attenuation values 
shifts downwards on the graph with decreasing film thickness. 

10 Figures 16(a) to 16(f) give contour plots of Re{S z ] related to the ss£ mode for 

films of thickness t = 20nm and widths w = 0.5\im and w = l^m, for three free-space 
wavelengths of operation: k 0 = 0.6, 0.8 and \.2\im. Comparing the contours shown in 
Figures 16(a) to 16(f), explains in part the frequency dependent behaviour plotted in 
Figures 15(a) and 15(b). Figures 16(a) to 16(f) show that the mode power contours 

15 spread out farther from the film as the wavelength increases, which means that the 
mode confinement to the metal region decreases, explaining in part the decrease in 
losses and the evolution of the mode towards the TEM mode of the background, as 
shown in Figures 15(a) and 15(b). This behaviour is more pronounced for the 
waveguide of width w = 0.5|im compared to the wider one of width w = \.Q\im. 

20 There are two mechanisms causing changes in the ss% mode as the frequency 

of operation varies. The first is geometrical dispersion, which changes the optical or 
apparent size of the film, and the second is material dispersion, which is modelled for 
the metal region using Equation (1). If no material dispersion is present, then the 
geometrical dispersion renders the film optically smaller as the free-space wavelength 

25 is increased (an effect similar to reducing t and w) so, in the case of the ssl mode, 
confinement to the film is reduced and the mode spreads out in all directions away 
from the latter. Now based on Equation (1), it is clear that the magnitude of the real 
part of the film's permittivity \Re{e r2 }\ varies approximately in a 1/co 2 or X\ fashion 
while the magnitude of its imaginary part |Iw{e f 2 }| varies approximately in a 1/w 3 

30 or A.J fashion. However, an increase in \Re{e r2 }\ reduces the penetration depth of the 
mode fields into the metal region and, combined with the geometrical dispersion, 
causes a net decrease in mode attenuation with increasing wavelength, even though the 
losses in the film increase in a k 0 fashion. 

Figure 15(b) shows that mode power attenuation values in the range 10 to 0. 1 

35 dB/cm are possible near communications wavelengths (k 0 ~ l.5\im) using structures of 
reasonable dimensions: w ~ l.0\xm and t ~ \5nm. Such values of attenuation are low 
enough to consider the ssl mode as bein 8 long-ranging, suggesting that these 
waveguides are practical for applications requiring propagation over short distances. 
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As shown in the previous section, even lower attenuation values are possible if the 
background permittivity is lowered. From Figures 16(e) and 16(f), (case X 0 = 1.2\im, 
which is near communications wavelengths), it is apparent that the mode power 
confinement is within one free- space wavelength of the film, which should be tight 
5 enough to keep the mode bound to the structure if a reasonable quality metal film of 
the right geometry can be constructed. 

VII. Mode Characteristics and Evolution With Film Thickness: Small Asymmetry 

A. Mode Solutions for a Metal Film Slab Waveguide 

10 Effects on waveguiding characteristics of using an asymmetric waveguide 

structure will now be discussed, beginning with the reproduction of results for an 
infinitely wide asymmetric metal film waveguide (similar to that shown in Figure 17(a) 
but with w = °°), taken from the standard work on such structures [6]. In order to 
remain consistent with their results, the optical free-space wavelength of excitation is 

15 set to X 0 = 0.633 (im and the value they used for the relative permittivity of the silver 
film at this wavelength is used here: e r2 = -19 - j0.53. The relative permittivities of 
the bottom and top dielectric regions are set to e T l = 4 (n, = 2) and e r 3 = 3.61 (n 3 = 
1.9); these values create a structure having a small asymmetry with respect to the 
horizontal dimension. 

20 The dispersion curves of the s b and a b modes supported by the infinitely wide 

structure were computed using the MoL and the results are shown in Figures 18(a) and 
18(b). From these figures, it is seen that the propagation constant of the a b mode tends 
towards that of the plasmon-polariton mode supported by the bottom interface, given 
•by Equations (14) and (15), as the thickness of the film increases. It is also noted that 

25 this mode does not exhibit a cutoff thickness, while it is clear that the s b mode has one 
near t = 18 nm. The propagation constant of the s b mode is seen to tend towards the 
value of a plasmon-polariton mode supported by the top interface as the thickness 
increases. These results are in perfect agreement with those reported in [6]. 

30 B. Modes Supported by a Metal Film of Width w = 1 pm 

The study proceeds with the analysis of the structure shown in Figure 17(a) for 
the case w = 1 [im. The material parameters and free-space wavelength that were used 
in the previous case w == «> were also used here. The dispersion curves for the first 
seven modes were computed using the MoL and the results are shown in Figures 18(a) 

35 and 18(b). 

In this asymmetric structure, true field symmetry exists only with respect to the 
y axis. With respect to the horizontal dimension, the modes have a symmetric-like or 
asymmetric-like field distribution with field localization along either the bottom or top 
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metal-dielectric interface. The modes that have a symmetric-like distribution with 
respect to the horizontal dimension are localized along the metal-dielectric interface 
with the lowest dielectric constant, while modes that have an asymmetric-like 
distribution with respect to this axis are localized along the metal-dielectric interface 
5 with the highest dielectric constant. This behaviour is consistent with that observed for 
asymmetric metal slab waveguides. 

The mode nomenclature adopted for symmetric structures can be used without 
ambiguity to describe the modes supported by asymmetric structures as long as the 
modes are identified when the metal film is fairly thick, before significant coupling 

10 begins to occur through the metal film, and while the origin of the mode can be 

identified unambiguously. As the metal film thickness decreases, the modes (and their 
fields) can evolve and change considerably more in an asymmetric structure compared 
to a symmetric one. The number of extrema in the main transverse electric field 
component of the mode is counted along the lateral dimension at the interface where 

15 the fields are localized. This number is then used in the mode nomenclature. 

It was observed in Section III that the modes supported by symmetric structures 
are in fact supermodes created from a coupling of "edge" and "corner" modes 
supported by each metal-dielectric interface defining the structure. As the thickness and 
width of the metal decrease, the coupling between these interface modes intensifies 

20 leading to dispersion and possibly evolution of the supermode. In asymmetric 

structures, the bound modes are also supermodes created in a similar manner, except 
that dissimilar interface modes may couple to each other to create the supermode. For 
instance, a mode having one field extremum along the top interface (along the top edge 
bounded by the corners) may couple with a mode having three extrema along the 

25 bottom interface. The main selection criterion determining which interface modes will 
couple to create the supermode is similarity in the value of their propagation constants. 
For all modes supported by an asymmetric structure, an apparent symmetry or 
asymmetry with respect to the horizontal dimension can still be observed in the corner 
modes. 

30 The sa^ aa l> ss l and as l modes are the fundamental modes supported by 

the structure. The sal and aal modes are comprised of coupled corner modes, • 
resembling the corresponding modes in a symmetric structure, except that the fields are 
localized near the substrate. These two modes do not change in character as the 
thickness of the film decreases. A narrowing of the metal film would eventually break 

35 the degeneracy observed in Figures 18(a) and 18(b). 

For a sufficiently large thickness (about 100 nm, for the present structure), 
the ssl and asl modes are comprised of coupled corner modes much like the 
corresponding modes in a symmetric structure except that the fields are localized near 
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the superstrate. As the thickness of the metal film decreases, both of these modes 
begin to evolve, changing completely in character for very thin films. Figures 19(a) to 
19(d) show the evolution of the E y field component related to the ss£ mode as the 
thickness of the film ranges from 100 nm (Figure 19(a)) to 40 nm (Figure 19(d)). It is 
5 clearly seen that the mode evolves from a symmetric-like mode having fields localized 
near the superstrate to an asymmetric-like mode having fields localized along the 
substrate-metal interface. A similar evolution is observed for the as% mode. This 
change in character is also apparent in their dispersion curves: they follow the general 
behaviour of a symmetric-like mode for large thicknesses but then slowly change to 

10 follow the behaviour of an asymmetric-like mode as the thickness decreases. Since the 
substrate dielectric constant is larger than the superstrate dielectric constant, the mode 
is v v pulled" from a symmetric-like mode to an asymmetric-like mode (having field 
localization at the substrate-metal interface) as the metal film becomes thinner. 
Figures 20(a) to 20(d) show the E y field component related to 

15 the ssl and sal m °des for two film thicknesses. From these Figures it is noted that 
the top and bottom edge modes comprising a supermode are different from each other. 
In Figure 20(a), for instance, it is seen that the bottom edge mode has three extrema 
and is of higher order than the top edge mode which has one extremum. A similar 
observation holds for Figure 20(c), where it can be seen that the bottom edge mode has 

20 one extremum while the top one has none. 1 In this structure, the substrate has a higher 
dielectric constant than the superstrate so the phase constant of a particular substrate- 
metal interface mode will be higher than the phase constant of the same mode at the 
metal-superstrate interface. Since a supermode is created from a coupling of edge 
modes having similar propagation constants, it should be expected that, in an 

25 asymmetric structure, different edge modes may couple to create a supermode. In 
general, higher-order modes have smaller values of phase constant compared to lower- 
order modes, so in structures having e 3 < € l5 all supermodes are comprised of a 
bottom edge mode of the same order or higher than the top edge mode, as shown in 
Figures 20(a) to 20(d). If e 3 > e b then the opposite statement is true. 

30 A careful inspection of the fields associated with the. ssl, sal and aal modes 

reveals that, as the thickness of the film decreases, the mode fields may evolve in a 
smooth manner similar to that shown in Figures 19(a) to 19(d), but, in addition, a 
change or "switch" of the constituent edge modes may also occur. For instance, from 
Figure 20(c), the sal mode is seen to comprise a substrate-metal interface mode 

35 having one extremum for a film thickness of 100 nm, while for a thickness of 60 nm 
the substrate-metal interface mode has three extrema, as shown in Figure 20(d). Since 
higher-order modes have in general lower phase constants than lower-order modes, this 
change in edge modes causes a reduction in the phase constant of the sa b mode in the 
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neighbourhood of 60 nm, as shown in Figure 18(a). Another change occurs near 40 nm 
as the corner modes switch from being symmetric-like (as in Figures 20(c) and 20(d)) 
to being asymmetric-like with respect to the horizontal dimension. This change is again 
reflected in the dispersion curve of the sal mode as its phase constant is seen to 
5 increase with a further decrease in thickness. In general, the changes in the edge and 
corner modes are consistent with the directions taken by the dispersion curves as the 
film thickness decreases, thus explaining the oscillations in the curves seen in Figures 
18(a) and 18(b). 

The only potentially long-ranging mode supported by this structure at the 
10 wavelength of analysis is the ssl mode. As shown in Figures 18(a) and 18(b), the 
mode has a cutoff thickness near t = 22 nm and although the attenuation drops quickly 
near this thickness, it should be remembered that the field confinement does so as well. 
Furthermore, the spatial distribution of the main transverse field component related to 
this mode evolves with decreasing thickness in the manner shown in Figure 20(a) and 
15 20(b), such that near cutoff the spatial distribution has strong extrema along the top and 
bottom edges. These extrema render the mode less excitable using an end-fire 
technique, so coupling losses would be higher compared to the fundamental symmetric 
mode in symmetric waveguides. Also, the fact that the mode would be operated near 
its cutoff thickness implies that very tight tolerances are required in the fabrication of 
20 structures. Nevertheless, it should be possible to observe propagation of this mode in a 
suitable structure using an end-fire experiment. 

VIII. Mode Characteristics and Evolution With Film Thickness: Large Asymmetry 

A. Mode Solutions for a Metal Film Slab Waveguide 

25 The study proceeds with the analysis of structures having a large difference in 

the dielectric constants of the substrate and superstate. With respect to Figure 17(a), 
the relative permittivities of the substrate and superstate are set to e r , = 4 (n x == 2) 
and e f 3 = 2.25 (n 3 = 1.5), respectively, the width of the metal film is set to w == <*>, 
and the dielectric constant of the metal region and the wavelength of analysis are set to 

30 the same values as in Section III. The dispersion curves of the s b and a b modes 

supported by this structure can be seen in Figures 21(a) and 21(b).. Comparing with 
Figures 18(a) and 18(b), it is observed that the s b mode has a larger cutoff thickness in 
a structure having a large asymmetry than in a structure having similar substrate and 
superstate dielectric constants. The results shown were computed using the MoL and 

35 are in perfect agreement with those reported in [6]. 



B. Modes Supported by a Metal Film of Width w == 7 ym 
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The structure shown in Figure 17(a) was analyzed using the MoL for w = 1 [im 
and for the same material parameters and free-space wavelength as those given above 
for w = °o. The dispersion curves of the first six modes supported by the structure are 
shown in Figures 21(a) and 21(b). 
5 An inspection of the mode fields related to the sal an ^ m l m °des reveals 

that these modes are again comprised of coupled corner modes with fields localized at 
the substrate-metal interface. The modes do not change in character as the thickness of 
the film decreases and a narrowing of the metal film would eventually break the 
degeneracy observed in Figures 21(a) and 21(b). 

10 The spatial distribution of the Ey field component related to the, ss£, asl, sal . 

and aal modes is given in Figures 22(a) to 22(d). It is noted from this figure that in 
all cases the metal-superstrate interface modes are similar: they have fields with no 
extrema along the interface but rather that are localized near the corners and have 
either a symmetric or asymmetric distribution with respect to the y axis. These corner 

15 modes are in fact the lowest order modes supported by the metal-superstrate interface; 
they have the largest value of phase constant and thus are most likely to couple with 
edge modes supported by the substrate-metal interface to form a supermode. From 
Figures 22(a) and 22(b) it is observed that the substrate-metal interface modes 
comprising the ssl and asl modes are of very high order. This is expected since the 

20 substrate dielectric constant is significantly higher than the superstrate dielectric 
constant and higher order modes have lower values of phase constant. 
The ssl and asl m °des shown in Figures 22(a) and 22(b) indeed have fields that are 
localized along the metal-superstrate interface, while the sal anc * m b modes shown 
in Figures 22(c) and 22(d) have fields that are localized along the substrate-metal 

25 interface. 

One effect, caused by increasing the difference between the substrate and 
superstrate dielectric constants, is that the difference between the orders of the top and 
bottom edge modes comprising a supermode can increase. This effect can be observed 
by comparing Figure 19(a) with Figure 22(a). In the former, there is no difference 

30 between the orders of the top and bottom edge modes, while in the latter the difference 
in the orders is 5. Another effect is that the degree of field localization increases near 
the interface between the metal and the dielectric of higher permittivity, for all modes 
that are asymmetric-like with respect to the horizontal dimension. This effect can be 
seen by comparing the fields related to the sal mode shown in Figures 22(c) and 

35 20(c). A comparison of the fields related to the sal and aal modes reveals that this 
effect is present in these modes as well. 

From the dispersion curves shown in Figure 21(a), it is apparent that the 
normalized phase constant of all modes converge with increasing film thickness to 
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normalized phase constants in the neighbourhood of those supported by plasmon- 
polariton waves localized along the associated isolated edge. The normalized phase 
constants of modes having fields localized at the substrate-metal interface, converge 
with increasing film thickness to normalized phase constants in the neighbourhood of 
5 that related to the a b mode, while the normalized phase constants of modes having 
fields localized along the metal-superstrate interface converge to values near that of the 
s b mode. This behaviour is present though less apparent, in structures where the 
asymmetry is smaller, such as the one analyzed in Section VII. 

Comparing Figures 18(a) and 18(b) withj Figures 21(a) and 21(b), it is noted 

10 that the dispersion curves of the modes are much smoother when the difference in the 
substrate and superstrate dielectric constants is large. This is due to the fact that the 
edge modes comprising the supermodes are less likely to change or switch as they do in 
a structure having similar substrate and superstrate dielectric constants. Thus modes 
that start out being symmetric-like with respect to the horizontal dimension remain so 

15 as the thickness of the film decreases. The cutoff thickness of the symmetric-like 
modes also increases as the difference between the substrate and superstrate dielectric 
constants increases. 

It is apparent that introducing a large asymmetry can hamper the ability of the 
structure to support useful long-ranging modes. Any mode that is long-ranging would 
20 likely have fields with numerous extrema along the width of the interface between the 
metal film and the dielectric of higher permittivity, as shown in Figures 22(a) and 
22(b). 

IX. Mode Dispersion with Film Width: Small Asymmetry 
25 An asymmetric structure comprising the same dielectrics as the structures 

studied in Section VII, but having a metal film of width w = 0.5 |im, was analyzed at 
the same free-space wavelength in order to determine the impact of a narrowing film 
width on the modes supported. The structure was analyzed using the MoL and Figures 
23(a) and 23(b) give the dispersion curves obtained for the first few modes supported. 
30 Comparing Figures 23(a) and 23(b) with Figures 18(a) and 18(b). reveals that 

reducing the width of the film does not cause major changes in the behaviour of the 
fundamental modes, but does have a major impact on the higher order modes. It is 
noted that reducing the film width increases the cutoff thickness of the ss£ mode. This 
higher order mode is symmetric-like with respect to the horizontal dimension, and the 
35 cutoff thickness of the symmetric-like modes in general increases as the width of the 
film decreases due to a reduction in field confinement to the metal film. 
The aa b mode was sought but not found for this film width. 
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It is also noted by comparing Figures 23(a) and 23(b) with Figures 18(a) and 
18(b) that the sal m( >de evolves quite differently depending on the width of the film. 
For a film width of w = 1 jim, the mode follows the general behaviour of an 
asymmetric-like mode whereas, for a film width of w = 0.5 fim, the mode evolves as 
5 a symmetric-like mode, and has a cutoff thickness near t = 27 nm. When the film is 
wide, it becomes possible for numerous higher order edge modes (having similar 
values of phase constant) to be supported by the substrate-metal or metal- superstate 
interfaces, so edge modes comprising a supermode are likely to change or switch as the 
thickness of the film is reduced, as shown in Figures 20(c) and 20(d). For a narrow 

10 metal film, some of the higher order edge modes may be cutoff, thus rendering 
changes in edge modes impossible. In such a case, the supermode may be forced to 
evolve in a smooth manner with decreasing film thickness. A close inspection of the 
mode fields related to the sal m °de for a film width of w = 0.5 |im reveals that there 
are no changes to the edge modes as the thickness decreases; rather the mode evolves 

15 smoothly from its field distribution at a large thickness (similar to that shown in Figure 
20(c)) to having a symmetric-like distribution with only one extremum along the top 
and bottom edges of the film. A change in behaviour due to a change in the width of 
the metal film was observed only for the sal m °de in this study, but such changes are 
in general not limited to this mode. 

20 The sal and ssl modes could be made to propagate over useful distances in 

this structure, if they are excited near their cutoff thicknesses. However, the difficulties 
outlined in Section VII B regarding the excitation of modes near cutoff also hold here. 

X. Evolution of the ss% and sal Modes With Structure Asymmetry 

25 The ss% and sal modes are of practical interest. The ss* mode is the main 

long-ranging mode supported by symmetric finite-width metal film structures, and, as 
demonstrated in the previous section, the sal mode can be the main long-ranging 
mode supported by asymmetric finite-width structures. In metal films of the right 
thickness, they are also the modes that are the most suitable to excitation in an end-fire 

30 arrangement. 

Structures comprising a substrate dielectric having = 2, of a metal film of 
width w = 0.5 |im, and of various superstate dielectrics having n 3 = 2, 1.99, 1.95 
and 1.9 were analyzed at the same free-space wavelength as in Section VII. The 
equivalent permittivity of the metal film was also set to the same value as in Section 
35 VII. These analyses were performed in order to investigate the effects on the 

propagation characteristics of the ss£ and sal mo d es caused by a slight decrease in 
the superstrate permittivity relative to the substrate permittivity. Figures 24(a) and 
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24(b) show the dispersion curves with film thickness, obtained for these modes in the 
four structures of interest. 

As seen in Figure 24(a) and its inset, the dispersion curves of the modes 
intersect at a certain film thickness only for the symmetric case (n 3 = n^. As soon as 
5 some degree of asymmetry exists, the curves no longer intersect, though they may 
come quite close to each other if the asymmetry is small, as seen in the case of n 3 = 
1.99. As the degree of asymmetry increases, the separation between the curves 
increases. 

The evolution with film thickness of the sal m °de is shown in Figures 25(a) 
10 to 25(d) for the case n 3 = 1.99 and for thicknesses about t = 59 nm (near the 

maximum in its phase dispersion curve). The evolution of this mode for the cases n 3 = 
1.95 and 1.9 is similar to that shown. The evolution with film thickness of 
the ss£ mode is similar in these structures to the evolution shown in Figures 19(a) to 
19(d) for the case w = 1 nm and n 3 = 1.9. Comparing Figures 25(a) to 25(d) and 
15 Figures 19(a) to 19(d), reveals that the modes v v swap" character near t = 59 nm. For 
film thicknesses sufficiently above this value, the modes exhibit their defining character 
as shown in Figures 19(a) and 25(a), but for film thicknesses below it, each mode 
exhibits the other's character, as shown in Figures 19(d) and 25(d). This character 
swap is present for the three cases of asymmetry considered here (n 3 =1.99, 1.95 and 
20 1.9) and explains the behaviour of the dispersion curves shown in Figures 24(a) and 
24(b). 

From Figures 24(a) and 24(b), it is noted that a cutoff thickness exists for the 
long-ranging mode as soon as an asymmetry is present in the structure. It is also 
observed that the cutoff thickness increases with increasing asymmetry. In the case of 

25 n 3 = 1.99, the cutoff thickness of the mode is near t = 12 nm, while for n 3 = 1.9 the 
cutoff thickness is near t = 27 nm. As the width of the metal film w increases, the 
cutoff thickness of the sal m °de decreases as long as the mode remains long-ranging 
(recall that the character of this mode may also change with film width such that its 
behaviour is similar to the a 5 mode in the corresponding slab structure, as shown in 

30 Figures 18(a) and 18(b)). Also, it is clear from Figures 23(a) and 23(b) that the cutoff 
thickness of the sal mode is greater than the cutoff thickness of the s b mode 
supported by the corresponding slab structure. These results imply that the long- 
ranging mode supported by a thin narrow metal film is more sensitive to differences in 
the superstate and substrate permittivities than the s b mode supported by the 

35 corresponding slab structure. This is reasonable in light of the fact that, in finite-width 
structures, the mode fields tunnel through the metal as in slab structures, but, in 
addition, the fields also wrap around the metal film. 
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Figure 24(b) shows that near cutoff, the attenuation of the sal mod e supported 
by an asymmetric structure drops much more rapidly than the attenuation of 
the ssl mode supported by a symmetric structure. Thus, a means for range 
extension, similar to that observed in asymmetric slab structures [7], exists for metal 
5 films of finite width, though the difficulties related to the excitation of a mode near its 
cutoff thickness, as described in Section VII B, also apply here. 

Figures 26(a) to 26(d) show contour plots of 91 {S z } associated with the long- 
ranging modes for the four cases of superstate permittivity considered. S z is the z- 
directed component of the Poynting vector and its spatial distribution is computed from 
10 the spatial distribution of the mode fields using: 

S z = (E x H; - E y H x yi (6) 

where H* y denotes the complex conjugate of H x y . Figure 26(a) shows the contour 
plot associated with the ss b mode supported by a symmetric structure (n 3 = n, = 2) 

15 of thickness t — 20 nm. Figures 26(b), (c) and (d) show contours associated with 
the sal mo <te f° r the three cases of structure asymmetry considered. The contour 
plots shown in Figures 26(b), (c) and (d) are computed for film thicknesses slightly 
above cutoff, representative of the thicknesses that would be used to observe these 
long-ranging modes experimentally. From those figures, it is noted that the contour 

20 plots become increasingly distorted and the fields increasingly localized at the metal- 
superstate interface as the degree of asymmetry in the structure increases. It is also 
apparent by comparing Figures 26(a) and 26(d), that in an end-fire experiment, less 
power should be coupled into the sal mode supported by the asymmetric structure 
with n 3 = 1.9, compared to the ssl mode supported by the symmetric structure. End- 

25 fire coupling losses seem to increase with increasing structure asymmetry. 

The high sensitivity of the long-ranging mode supported by thin metal films of 
finite width, to structure asymmetry, is potentially useful. A small induced asymmetry 
(created via an electro-optic effect present in the dielectrics say) can evidently effect a 
large change in the propagation characteristics of the long-ranging mode. From 

30 Figures 24(a) and 24(b), it is apparent that a difference between the substrate and 
superstate refractive indices as small as n x - n 3 = An = 0.01 is sufficient to create an 
asymmetric structure where the long-ranging mode has a cutoff thickness of about t = 
12 nm. From Figure 24(a), a slightly larger difference of An = 0.05 changes the 
normalized phase constant of the long-ranging mode by A(p/p 0 ) ~ 0.025 for a metal 

35 film thickness of t = 20 nm. Both of these effects are potentially useful. 

Asymmetric structures having superstrate dielectric constants that are slightly 
greater than that of the substrate were also analyzed. The substrate dielectric constant 
was set to n, == 2 and superstrate dielectrics having n 3 = 2.01, 2.05 and 2.1 were 
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considered for the same metal, film width and operating wavelength. The results are 
similar to those presented in Figures 24 through 26 and the cut-off thicknesses are near 
those shown in Figure 24(b). Though the results are not identical, there is no major 
difference between the behaviour of the ss£ and sal modes 
5 whether e t > e 2 or z x < e 3 as long as the permittivities are similar. 

XI. Conclusion 

The purely bound optical modes supported by thin lossy metal films of finite 
width, embedded in an "infinite" homogeneous dielectric have been characterized and 

10 described. The modes supported by these symmetric structures are divided into four 
families depending on the symmetry of their mode fields and none of the modes are 
TM in nature (as they are in the metal film slab waveguide). In addition to the four 
fundamental modes that exist, numerous higher order modes are supported as well. A 
proposed mode nomenclature suitable for identifying them has been discussed. The 

15 dispersion of the modes with film thickness has been assessed, and the behaviour in 
general terms found to be consistent with that of the purely bound modes supported by 
the metal film slab waveguide. In addition, it has been found that one of the 
fundamental modes and some higher order modes have cut-off thicknesses. Mode 
dispersion with film width has also been investigated and it has been determined that 

20 the higher order modes have a cut-off width, below which they are no longer 

propagated. The effect of varying the background permittivity on the modes has been 
investigated as well, and the general behaviour found to be consistent with that of the 
modes supported by a metal film slab waveguide. In addition it was determined that 
the cut-off width of the higher order modes decreases with decreasing background 

25 permittivity and that all cut-off thicknesses are increased. 

One of the fundamental modes supported by the symmetric structures, the 
ssl mode, exhibits very interesting characteristics and is potentially quite useful. 
This mode evolves with decreasing film thickness towards the TEM wave supported by 
the background, (an evolution similar to that exhibited by the s b mode in metal film 

30 slab waveguides), its losses and phase constant tending asymptotically towards those of 
the TEM wave. In addition, it has been found that decreasing the film width can 
reduce the losses well below those of the s b mode supported by the corresponding metal 
film slab waveguide. Reducing the background permittivity further reduces the losses. 
However, a reduction in losses is always accompanied by a reduction in field 

35 confinement to the waveguide core, which means that both of these parameters must be 
traded-off one against the other. Furthermore, carefully selecting the film's thickness 
and width can make the ssl mode the only long-ranging mode supported. It has also 
been demonstrated that mode power attenuation values in the range of 10 to 0. 1 dB/cm 
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are achievable at optical communications wavelengths, with even lower values 
possible. Finally, evolved into its most useful form, the ss£ mode has a field 
distribution that renders it excitable using end-fire techniques. 

The existence of the ss£ mode in a symmetric structure, as well as its 
5 interesting characteristics, makes the finite-width metal film waveguide attractive for 
applications requiring short propagation distances. The waveguide offers two- 
dimensional field confinement in the transverse plane, rendering it useful as the basis 
of an integrated optics technology. Interconnects, power splitters, power couplers and 
interferometers could be built using the guides. Finally, the structures are quite simple 

10 and so should be inexpensive to fabricate. 

The long-ranging modes supported by asymmetric structures of finite width 
have a rapidly diminishing attenuation near their cutoff thickness (like asymmetric slab 
structures). The rate of decrease of the attenuation with decreasing thickness near 
cutoff is greater than the rate related to the ss£ mode in symmetric structures. 

15 However field confinement also diminishes rapidly near cutoff, implying that the 
structures ought to be fabricated to very tight tolerances and that all metal-dielectric 
interfaces should be of the highest quality. It has also been found that decreasing the 
width of the film increases the cutoff thickness of the main long-ranging mode. Below 
this cutoff thickness, no purely bound long-ranging mode exists. The long-ranging 

20 modes supported by metal films of finite-width are thus more sensitive to the 
asymmetry in the structure as compared to the s b mode supported by similar slab 
waveguides. This is a potentially useful result in that a small induced change in 
substrate or superstrate refractive index can have a greater impact on the long-ranging 
mode supported by a finite-width structure as compared to a similar slab waveguide. 

25 Parts of the foregoing theoretical discussion have been published by the 

inventor in references [13] [20], [44] and [45]. 

SPECIFIC EMBODIMENTS AND EXAMPLES OF APPLICATION 

Examples of practical waveguide structures, and integrated optics devices which 
30 can be implemented using the invention, will now be described with reference also to 
Figures 27 to 42, 51 and 52. Unless otherwise stated, where a waveguide structure is 
shown, it will have a general construction similar to that shown in Figures 1(a) and 
1(b) or that shown in Figures 17(a) and 17(b). 

Waveguides, structures and devices disclosed herein operate with radiation: 
35 • having a wavelength such that a plasmon-polariton wave is supported; 

• at optical wavelengths; 

• at optical communications wavelengths; 

• at wavelengths in the range of 800 nm to 2 [im; 
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• at wavelengths near 1550 nm; 

• at wavelengths near 1310 nm; 

• at wavelengths near 850 nm; 

• at wavelengths near 980 nm. 

5 

It should be appreciated that references to wavelength should be interpreted as 
meaning the centre free space wavelength of the spectrum associated with the input 
radiation. 

The waveguide structure 100 shown in Figures 1(a) and 1(b) comprises a strip 

10 of finite thickness t and width w of a first material having a high free (or almost free) 
charge carrier density, surrounded by a second material which has a very low free (or 
almost free) charge carrier density. The strip material can be a metal or a highly doped 
semiconductor and the background material can be a dielectric, or an undoped or low 
doped semiconductor. The thickness t and the width w of the strip are selected such 

15 that the waveguide supports the main long-ranging plasmon-polariton mode of interest, 
the ss b ° mode, at the free-space wavelength of interest. 

The waveguides will propagate plasmon-polariton waves if the strip has a width 
in the range from about 0. 1 \im to about 12 \im and a thickness in the range from about 
5 nm to about 100 nm, particular dimensions depending on the index of refraction of 

20 the surrounding material, the strip material and the free- space wavelength of operation. 
Figures 51 and 52 give the mode power attenuation for example waveguides 
constructed from strips of gold (Au) and aluminium (Al), respectively, each embedded 
in silicon dioxide (SiO a ), for various widths and thicknesses of the metal strip. The 
analyses were carried out at an optical free space wavelength of 1550 nm. The curves 

25 show that very low attenuation values can be obtained with metal strips of practical 
dimensions. Generally, the attenuation using the gold strip is about one half of that 
obtained with the aluminium strip having similar dimensions. Suitable ranges for the 
dimensions of the strip, in these particular waveguide examples, can be deduced from 
Figures 51 and 52. A suitable range for the thickness of the strip is from about 5 nm 

30 to about 50 nm, and a suitable range for the width of the strip is from about 0.5 jim to 
about 12 \irn. For either of these metals, particularly good dimensions are a thickness 
of about 20 nm and a width of about 4 (am; such a strip yields a low loss waveguide 
having a mode size of about 10 jim. 

For a waveguide structure having particular materials for the strip and 

35 surrounding material, operable at a particular wavelength, there will be a range of 
thickness and width for the strip which will allow the plasmon-polariton wave to be 
supported. For metallic strips surrounded by material such as silicon dioxide, silicon 
oxynitride or silicon nitride, having a refractive index in the range from about 1.4 to 
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about 2.0, the range of dimensions for the metallic strip is the width in the range from 
about 0.5 |im to about 12 \im and the thickness in the range from about 5 nm to about 
50 nm. Such waveguide structures support effectively the propagation of a plasmon- 
polariton wave having a wavelength in the range from about 0.8 \xm to about 2 jim. 
5 For operation in the wavelength range from about 1.3 \im to about 1.7 \im, good 
dimensions for the metallic strip are a width in the range from about 0.7 \im to about 8 
\xm and a thickness in the range from about 15 nm to about 25 nm. A particularly good 
choice for the dimensions of the metallic strip is a width of about 4 |im and a thickness 
of about 20 nm for operation at a wavelength near 1.55 \im. Other materials having an 

10 index of refraction approximately in the same range require strip dimensions in 
approximately the same ranges. 

For strips surrounded by materials such as lithium niobate or PLZT, having a 
refractive index in the range from about 2 to about 2.5, and required to support 
effectively the propagation of a plasmon-polariton wave having a free space wavelength 

15 in the range from about 0.8 |im to about 2 nm, the ranges of dimensions for the strip 
preferably include a width in the range from about 0. 15 \xiu to about 6 |im and a 
thickness in the range from about 5 nm to about 80 nm. For operation in the 
wavelength range from about 1.3 \im to about 1.7 \im, good dimensions for the strip 
are a width in the range from about 0.4 |im to about 2 \im and thickness in the range 

20 from about 15 nm to about 40 nm. For operation at a wavelength near 1.55 \im, the 
dimensions of the strip preferably are a width of about 1 jim and a thickness of about 
20 nm. Other surrounding materials having an index of refraction approximately in the 
same range require strip dimensions in approximately the same ranges. 

For strips surrounded by materials such as silicon, gallium arsenide, indium 

25 phosphide or other material having a refractive index in the range from about 2.5 to 
about 3.5, and required to support effectively the propagation of a plasmon-polariton 
wave having a free space wavelength in the range from about 0.8 fim to about 2 jim, 
the ranges of dimensions for the strip preferably include a width in the range from 
about 0. 1 |im to about 2 \im and thickness in the range from about 5 nm to about 30 

30 nm. For operation in the wavelength range from about 1.3 \im to about 1.7 ^m, good 
dimensions for the strip are a width in the range from about 0. 13 |im to about 0.5 |im 
and thickness in the range from about 10 nm to about 20 nm. For operation at a 
wavelength near 1.55 nm, the dimensions. of the strip preferably are a width of about 
0.25 \im and thickness of about 15 nm. Other surrounding materials having an index of 

35 refraction approximately in the same range require strip dimensions in approximately 
the same ranges. 
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Unless otherwise stated, when waveguide structure dimensions are mentioned 
from this point onward, they refer to the Au/Si0 2 material combination at an operating 
optical free- space wavelength of 1550 nm. 

The plasmon-polariton wave which propagates along a waveguide structure may 
5 be excited by an appropriate optical field incident at one of the ends of the waveguide, 
as in an end-fire configuration, and/or by different radiation coupling means. In order 
to achieve good coupling efficiency with the input and output means, the plasmon- 
polariton wave preferably has a transverse electric field component that overlaps well 
with that of the fundamental mode of propagation supported by the input and output 

10 means. Tapered input and output strip sections can be used to match the mode sizes of 
the plasmon-polariton waveguide and the waveguide used as the input and output 
means, in order to achieve high coupling efficiency . An example of an appropriate 
end-fire arrangement comprises a standard optical fibre butt-coupled to the input of the 
waveguide (the output of the waveguide can also be butt-coupled to a fibre). 

15 An example waveguide structure achieving high coupling efficiency with 

standard optical fibre at a free space wavelength near 1.55 [im comprises an Au strip 
surrounded by Si0 2 , the strip having a thickness of about 20 nm and a width in the 
range from about 4 |im to about 8 |im. Alternative means for exciting the waveguide 
include excitation at an intermediate position using, for example, a prism and the so- 

20 called attenuated total reflection method (ATR). 

Suitable materials for the surrounding material include (but are not limited to) 
glass, glasses, quartz, polymer, silicon dioxide, silicon nitride, silicon oxynitride, and 
undoped or very lightly doped GaAs, InP or Si. The material may comprise single 
crystal or partially crystalline material. The surrounding material may comprise 

25 electro-optic crystals or material such as lithium niobate (LiNb0 3 ), PLZT or electro- 
optic polymers. The surrounding material is not necessarily homogeneous. 

Suitable materials for the strip include (but are not limited to) metals, semi- 
metals, and highly n- or p-doped semiconductors such as GaAs, InP or Si. Suitable 
metals for the strip may comprise, a single metal, or a combination of metals, selected 

30 from the group Au, Ag, Cu, Al, Pt, Pd, Ti, Ni, Mo, and Cr, preferred metals being 
Au, Ag, Cu, and Al. Metal silicides such as CoSi 2 are particularly suitable for use with 
Si as the surrounding material. A single material or combination of materials which 
behave like metals, such as Indium Tin Oxide (ITO) can also be used. In the context of 
this specification, the term "metallic" substances embraces such material(s). The 

35 metallic strip is not necessarily homogeneous. 

Particularly suitable combinations of materials include Au for the strip and Si0 2 
for the surrounding material, or Au for the strip and LiNb03 for the surrounding 
material, or Al or CoSi 2 for the strip and Si for the surrounding material, or Au for the 
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strip and polymer for the surrounding material, or gold for the strip and PLZT for the 
surrounding material. These combinations are preferred because, in each case, the 
strip material exhibits low optical absorption and is chemically compatible with the 
surrounding material. Devices may be fabricated using known deposition techniques or 
5 wafer bonding and polishing for the cladding materials, and known lithographic and 
deposition techniques for the strip. 

The waveguide structures may comprise: a strip that is homogeneous and a 
surrounding material that is homogeneous, a strip that is homogeneous and a 
surrounding material that is inhomogeneous; a strip that is inhomogeneous and a 

10 surrounding material that is homogeneous; or a strip that is inhomogeneous and a 
surrounding material that is inhomogeneous. An inhomogeneous strip can be formed 
from a continuously variable material composition, or strips or laminae. An 
inhomogeneous surrounding material can be formed from a continuously variable 
material composition, or strips or laminae. 

15 Generally, a plasmon-polariton waveguide having a metallic strip of large 

aspect ratio supports substantially TM polarized light, ie.; the transverse electric field 
component of the plasmon-polariton wave, the ss b ° mode, is aligned substantially with 
the normal to the largest dimension in the waveguide cross-section. Thus, a straight 
waveguide 100 with the dimensions set out above is polarisation sensitive. The 

20 plasmon-polariton wave is highly linearly polarised in the vertical direction, i.e. 
perpendicular to the plane of the strip. Hence, it may serve as a polarisation filter, 
whereby substantially only a vertical polarised mode (aligned along the y-axis as 
defined in Figure 1(a)) of the incident light is guided. 

The length / shown in Figure 1(b) is arbitrary and will be selected to implement 

25 a desired interconnection. 

Figure 27 shows a transition waveguide section 102 having stepped sides which 
can be used to interconnect two sections of waveguide having different widths. The 
larger width can be used to more effectively couple the waveguide to input/output 
means, such as fibres. The reduced width helps to reduce the insertion loss of the 

30 waveguide. Preferred widths are about W 2 = 8 \im to couple to single mode fibre and 
about W, =4 |am for the waveguide width. Any symmetry of the structure shown can 
be used. 

Figure 28 shows an angled section 104 which can be used as an interconnect. 
Its dimensions, W l5 W 2 and 1 and the angles (J> t and (J) 2 , are adjusted for a particular 
35 interconnection application as needed. Usually the angles are kept small, in the range 
of about 1 to about 15 degrees and the input and output widths are usually similar, 
about 4 fim. Although the sides of the angled section 104 shown in Figure 28 are 
tapered, they could be parallel. It should also be appreciated that the angle of the 



42 

inclination could be reversed, i.e. the device could be symmetrical about the bottom 
right hand corner shown in Figure 28 or transposed about that axis if not symmetrical 
about it. 

Figure 29 shows a tapered waveguide section 106, which can be used to 
5 interconnect two waveguides of different widths. The length of the taper is usually 
adjusted such that the angles are small, preferably in the range from about 0.5 to about 
15 degrees. The taper angles at the two sides are not necessarily the same. Such a 
configuration might be used as an input port, perhaps as an alternative to the layout 
shown in Figure 27, or as part of another device, such as a power splitter. The larger 
10 width can be used to more effectively couple the waveguide to input/output means, 
such as fibres. The reduced width helps to reduce the insertion loss of the waveguide. 
Typical widths are about W 2 =8 (im to couple to standard single mode fibre and about 
W, =4 \xm for the waveguide width. Any symmetry of the structure shown can be 
used. 

15 Figure 30 illustrates an alternative transition waveguide section 130 which has 

curved sides, rather than straight as in the trapezoidal transition section disclosed in 
Figure 29. In Figure 30, the curved sides are shown as sections of circles of radius 
and R 2 , subtending angles <J>i and 4> 2 respectively, but it should be appreciated that 
various functions can be implemented, such as exponential or parabolic, such that the 

20 input/output reflections and the transition losses are minimised. 

Figure 31 shows a curved waveguide section 108 which can be used to redirect 
the plasmon-polariton wave. The angle $ of the bend can be any angle suitable for a 
particular implementation. 

In addition to the width and thickness of the strip, critical dimensions in the 

25 design of a bend include the bending radius R and offsets O, and 0 2 determining the 
positions of the input and output straight sections 100. Although the bend will work 
for many strip dimensions and bend designs, i.e., the structure 108 will convey the 
plasmon-polariton wave around the bend, there may be (i) radiation loss or leakage out 
of the bend, (ii) transition losses at the planes where the bent section meets the straight 

30 sections, and possibly also (Hi) reflections back towards the direction from which the 
wave came. Minimal radiation loss or leakage is obtained by selecting a radius of 
curvature for the strip that is greater than a threshold value, for which radiation is 
negligible. 

Reduced transition and reflection losses may be obtained by offsetting the input 
35 and output waveguides 100 radically outwards relative to the ends of the bend. The 
reason for this is that the straight waveguide sections 100 have an optical field 
extremum that peaks along the longitudinal centre line of the strip, and decays towards 
its edges. In the bend, the extremum of the optical field distribution shifts towards the 
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exterior of the bend. This results in increased radiation or leakage from the bend and 
increased transition and reflection losses due to a mismatch in the field distributions. 
Offsetting the input and output waveguides 100 towards the outside of the bend aligns 
the extrema of their optical fields more closely with that of the optical field in the bent 
5 section 108, which helps to reduce, even minimise, both transition and reflection 
losses. The smaller the radius R, the greater the radiation from the bend. Also the 
offsets O, and 0 2 are related to the radius R. Optimum values for these parameters 
would have to be determined according to the specific bending requirements. 

It should also be noted that it is not necessary to connect the input and output 

10 waveguides 100 directly to the curved section. As shown in Figure 31, it is possible to 
have a short spacing d, between the end of the input waveguide 100 and the adjacent 
end of the curved section 108. Generally speaking, that spacing d r should be 
minimised, even zero, and probably no more than a few optical wavelengths. A 
similar offset 0 2 and spacing d 2 could be provided between the bend 108 and the output 

15 straight waveguide 100. 

Although Figure 31 shows no gradual transition between the straight 
waveguides 100 at the input and output and the ends of the curved section 108, it is 
envisaged that, in practice, a more gradual offset could be provided so as to reduce 
edge effects at the corners. 

20 Estimates of radiation loss associated with 90° bends were obtained as a 

function of the bend's radius of curvature R for various strip widths and thickness 
waveguide bends comprised of a gold (Au) strip, and for two cases of surrounding 
material: silicon dioxide (Si0 2 ) and lithium niobate (LiNb03). The results clearly 
indicated the onset of significant radiation loss as the radius of curvature for a 

25 particular strip geometry was reduced. It was inferred from the results that good 
bending radii are in the range from about 100 \im to about 10 cm for a strip width in 
the range from about 0. 1 |xm to about 12 |im and a strip thickness in the range from 
about 5 nm to about 100 nm. 

Acceptable 90 degree bends, comprising strips surrounded by materials such as 

30 silicon dioxide, silicon oxynitride or silicon nitride, having a refractive index in the 
range from about 1 .4 to about 2.0, the range of dimensions for the metallic strip being 
a width in the range from about 0.5 |im to about 12 |im and a thickness in the range 
from about 5 nm to about 50 nm should have a radius of curvature in the range from 
about 100 [im to about 10 cm. Such bends support effectively the propagation of a 

35 plasmon-polariton wave having a wavelength in the range from about 0.8 |im to about 
2 |im. For operation in the wavelength range from about 1.3 |im to about 1.7 \im, 90 
degree bends have the following dimensions for the metallic strip: a width in the range 
from about 0.7 |im to about 8 [im, a thickness in the range from about 15 nm to about 
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25 nm and a radius of curvature in the range from about 1 mm to about 10 cm. For 
operation at a wavelength near 1 .55 nm, a particularly good choice for a 90 degree 
bend is a metallic strip having a width of about 6 |im, a thickness of about 20 nm and a 
radius of about 2 cm. Other materials having an index of refraction approximately in 
5 the same range require strip dimensions and radii of curvature in approximately the 
same ranges. 

For 90 degree bends comprising strips surrounded by materials such as lithium 
niobate or PLZT, having a refractive index in the range from about 2 to about 2.5, 
required to support effectively the propagation of a plasmon-polariton wave having a 

10 wavelength in the range from about 0.8 \xm to about 2 |im, the ranges of dimensions 
for the metallic strip include width in the range from about 0. 15 |im to about 6 nm, 
thickness in the range from about 5 nm to about 80 nm, and radius of curvature in the 
range from about 100 jam to about 10 cm. For operation in the wavelength range from 
about 1.3 \im to about 1.7 |im, 90 degree bends preferably have the following 

15 dimensions for the metallic strip: a width in the range from about 0.4 \xm to about 2 
nm, a thickness in the range from about 15 nm to about 40 nm and a ; radius of 
curvature in the range from about 1 mm to about 10 cm. A particularly good choice for 
a 90 degree bend is a metallic strip having a width of about 1.5 |im, a thickness of 
about 20 nm and a radius of about 4 cm for operation at a wavelength near 1 .55 \im. 

20 Other materials having an index of refraction approximately in the same range require 
strip dimensions and radii of curvature in approximately the same ranges. Similar 
ranges apply to semiconductors having a refractive index in the range from about 2.5 
to about 3.5, such as silicon, gallium arsenide and indium phosphide. 

The radius of curvature selected can be smaller than that used for a 90 degree 

25 bend if the desired bending angle is smaller than 90 degrees. 

Figure 32 shows a two-way power splitter 1 10 formed from a trapezoidal 
section 106 with a straight section 100 coupled to its narrower end 112 and two angled 
sections 104 coupled side-by-side to its wider end 1 14. The distances between the 
input waveguide 100 and the narrower end 112 of the tapered section 106, and between 

30 the output waveguides 104 and the wider end 1 14 of the tapered section 106, d u d 2 and 
d 3 , respectively, should be minimised. The angle between the output waveguides 104 
can be in the range of 0.5 to 10 degrees and their widths are usually similar. The 
offsets S t and S 2 between the output waveguides and the longitudinal centre line of the 
trapezoidal section 106 preferably are set to zero, but could be non-zero, if desired, 

35 and vary in size. Ideally, however, the output sections 104 should together be equal in 
width to the wider end 114. 

Offset Sj need not be equal to offset S 2 but it is preferable that both are set to 
zero. The widths of the output sections 104 can be adjusted to vary the ratio of the 
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output powers. The dimensions of the centre tapered section 106 are usually adjusted 
to minimise input and output reflections and radiation losses in the region between the 
output sections 104. 

It should also be noted that the centre tapered section 106 could have angles that 
5 vary according to application and need not be symmetrical. 

It is envisaged that the tapered section 106 could be replaced by a rectangular 
transition section having a width broader than the width of the input waveguide 100 so 
that the transition section favoured multimode propagation causing constructive and 
destructive interference patterns throughout its length. The length could be selected so 
10 that, at the output end of the rectangular transition section, the constructive portions of 
the interference pattern would be coupled into the different waveguides establishing, in 
effect, a 1-to-N power split. Such a splitter then would be termed a multimode 
interferometer-based power divider. 

It should be appreciated that the device shown in Figure 32 could also be used 
15 as a combiner. In this usage, the light would be injected into the waveguide sections 
104 and combined by the tapered centre section 106 to form the output wave which 
would emerge from the straight waveguide section 100. 

In either the Y splitter or the interferometer-based power divider, the number of 
arms or limbs 104 at the output could be far more than the two that are shown in 
20 Figure 32. 

It is also feasible to have a plurality of input waveguides. This would enable an 
NxN divider to be constructed. The dimensions of the transition section 106 then 
would be controlled according to the type of splitting/combining required. 

As shown in Figure 33(a), an angled waveguide section 104 may be used to 

25 form an intersection between two straight waveguide sections 100, with the dimensions 
adjusted for the particular application. It should be noted that, as shown in Figure 
33(a), the two straight sections 100 are offset laterally away from each other by the 
distances O, and 0 2 , respectively, which would be selected to optimise the couplings 
by reducing radiation, transition and reflection losses. The angle of the trapezoidal 

30 section 104 relative to the sections 100 will be a factor in determining the best values 
for the offsets 0 1 and 0 2 . The sections 100 and 104 need not be connected directly 
together, but could be spaced by the distances d! and d 2 and/or coupled by a suitable 
transition piece that would make the junction more gradual (i.e., the change of 
direction would be more gradual). 

35 The embodiments of Figures 31 and 33(a) illustrate a general principle of 

aligning optical fields, conveniently by offsets, wherever there is a transition or change 
of direction of the optical wave and an inclination relative to the original path, which 
can cause radiation, transition and reflection losses if fields are misaligned. Such 
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offsets could be applied whether the direction-changing sections were straight or 
curved. 

As shown in Figure 33(b), an alternative embodiment to that shown in Figure 
33(a) is obtained by replacing the angled section 104 by a pair of concatenated bend 
5 sections 108 arranged to form an S bend (with offsets) in order to provide the 
interconnection between the sections 100. The bend sections could be replaced with 
bend sections having a continuously variable radius of curvature implementing a bend 
profile other than circular. 

As illustrated in Figure 34(a), a power divider 1 16 can also be implemented 

10 using a pair of concatenated bend sections 108 instead of each of the angled sections 
104 in the splitter 1 10 shown in Figure 32. As shown in Figure 34(a), in each pair, 
the bend section nearest to the wider end 114 of the tapered section 106 bends outwards 
from the longitudinal centre line of the tapered section 106 while the other section 
bends oppositely so that they form an "S" bend. Also, the bend sections in each pair 

15 are offset by distance O x or 0 2 one relative to the other for the reasons discussed with 
respect to the bend 108 shown in Figure 31. Other observations made regarding the 
power divider and the bend section disclosed in Figures 31 and 32 respectively, also 
apply to this case. 

Figure 34(b) shows a modification to the power divider of Figure 34(a) 

20 specifically the use of a bifurcated transition section 106' having a narrower end 
portion 115A and two curved sections 115B and 115C diverging away from the 
narrower end portion 115 A; in effect, the transition section comprising two mirrored 
and overlapping bend sections. The light emerges from respective distal ends of the 
curved sections 1 15B and 115C. Two additional bend sections 108 are shown 

25 connected to distal ends of the two curved sections 1 15B and 1 15C, respectively, in 
effect forming with the transition section 106* two mirrored and overlapping S-bends. 

Figure 35 illustrates a Mach-Zender interferometer 118 created by 
interconnecting two power splitters 1 10 as disclosed in Figure 32. Of course, either or 
both of them could be replaced by the power splitter 116 shown in Figure 34(a) or 

30 Figure 34 (b). Light injected into one of the ports, i.e. the straight section 100 of one 
power splitter 110/116, is split into equal amplitude and phase components that travel 
along the angled arms 104 of the splitter, are coupled by straight sections 100 into the 
corresponding arms of the other splitter, and then are recombined to form the output 
wave. 

35 If the insertion phase along one or both arms of the device is modified, then 

destructive interference between the re-combined waves can be induced. This induced 
destructive interference is the basis of a device that can be used to modulate the 
intensity of an input optical wave. The lengths of the arms 100 are usually adjusted 
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such that the phase difference in the re-combined waves is 180 degrees for a particular 
relative change in insertion phase per unit length along the arms. The structure will 
thus be optically long if the mechanism used to modify the per unit length insertion 
phase is weak (or optically short if the mechanism is strong). 
5 Figure 36(a) illustrates a modulator 120 based on the Mach-Zender 118 

disclosed in Figure 35. As illustrated also in Figure 36(b), parallel plate electrodes 122 
and 124 are disposed above and below, respectively, each of the strips 100 which 
interconnects two angled sections 104, and spaced from it, by the dielectric material, at 
a distance large enough that optical coupling to the electrodes is negligible. The 

10 electrodes are connected in common to one terminal of a voltage source 126, and the 
intervening strip 100 is connected using a minimally invasive contact to the other 
terminal. Variation of the voltage V applied by source 126 effects the modulating 
action. According to the plasma model for the strip 100, a change in the carrier 
density of the latter (due to charging +2Q or -2Q) causes a change in its permittivity, 

15 which in turn causes a change in the insertion phase of the arm. (The change induced 
in the permittivity is described by the plasma model representing the guiding strip 100 
at the operating wavelength of interest. Such model is well known to those of ordinary 
skill in the art and so will not be described further herein. For more information the 
reader is directed to reference [21], for example.) This change is sufficient to induce 

20 destructive interference when the waves in both arms re-combine at the output 
combiner. 

Figure 36(c) illustrates an alternative connection arrangement in which the two 
plate electrodes 122 and 124 are connected to respective ones of the terminals of the 
voltage source 126. In this case, the dielectric material used as the background of the 
25 waveguide is electro-optic (LiNb0 3 , a polymer,...). In this instance, the applied 
voltage V effects a change in the permittivity of the background dielectric, thus 
changing the insertion phase along the arm. This change is sufficient to induce 
destructive interference when the waves in both arms re-combine at the output 
combiner. 

30 It will be noted that, in Figure 36(a), one voltage source supplies the voltage V! 

while the other supplies the voltage V 2 . W x and V 2 may or may not be equal. 

For both cases described above, it is possible to apply voltages in opposite 
polarity to both arms of the structure. This effects an increase in the insertion phase of 
one arm and a decrease in the insertion phase of the other arm of the Mach-Zender (or 
35 vice versa), thus reducing the magnitude of the voltage or the length of the structure 
required to achieve a desirable degree of destructive interference at the output. 

Also, it is possible to provide electrodes 122 and 124 and a source 126 for only 
one of the intervening strips 100 in order to provide the required interference. 
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It should be appreciated that other electrode structures could be used to apply 
the necessary voltages. For example, the electrodes 122 and 124 could be coplanar 
with the intervening strip 100, one on each side of it. By carefully laying out the 
electrodes as a microwave waveguide, a high frequency modulator (capable of 
5 modulation rates in excess of 10 Gbit/s) can be achieved. 

Figure 37 illustrates an alternative implementation of a Mach-Zender 128 which 
has the same set of waveguides as that shown in Figure 35 but which makes use of 
magnetic fields B applied to either or both of the middle straight section arms to induce 
a change in the permittivity tensor describing the strips. (The change induced in the 

10 tensor is described by the plasma model representing the guiding strip at the operating 
wavelength of interest. Such model is well known to those of ordinary skill in the art 
and so will not be described further herein. For more information the reader is 
directed to reference [21], for example.) The change induced in the permittivity tensor 
will induce a change in the insertion phase of either or both arms thus inducing a 

15 relative phase difference between the light passing in the arms and generating 
destructive interference when the waves re-combine at the output combiner. 
Modulating the magnetic field thus modulates the intensity of the light transmitted 
through the device. The magnetic field B can be made to originate from current- 
carrying wires or coils disposed around the arms 100 in such a manner as to create the 

20 magnetic field in the desired orientation and intensity in the optical waveguides. The 
magnetic field may have one or all of the orientations shown, B x , B y or B z or their 
opposites. The wires or coils could be fabricated using plated via holes and printed 
lines or other conductors in known manner. Alternatively, the field could be provided 
by an external source, such as a solenoid or toroid having poles on one or both sides of 

25 the strip. 

Figure 38 illustrates a periodic waveguide structure 132 comprising a series of 
unit cells 134, where each cell 134 comprises two rectangular waveguides 100 and 100' 
having different lengths /, and l 2 and widths w ( and w 2 , respectively. The dimensions 
of the waveguides in each unit cell 134, the spacing therebetween, the number of 

30 unit cells, and the spacings d 2 between cells are adjusted such that reflection occurs at a 
desired operating wavelength or over a desired operating bandwidth for an optical 
signal propagating along the grating axis, i.e. the longitudinal axis of the array of cells 
134. The period of the periodic structure, i.e. the length of each unit cell, /, + l 2 + d l 
+ d 2 , can be made optically long, such that a long-period periodic structure is 

35 obtained. The dimensions of the elements 100, 100' in each unit cell 134 can also be 
made to change along the direction of the periodic structure in order to implement a 
prescribed transfer function (like in a chirped periodic structure). 
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It should be noted that the waveguides 100, 100' in each cell need not be 
rectangular, but a variety of other shapes could be used. For example, Figure 39 
illustrates a portion, specifically two unit cells 138 only, of an alternative periodic 
structure 136 in which each unit cell 138 comprises two of the trapezoidal waveguide 
5 sections 106, 106' like that described with reference to Figure 30, with their wider 
edges opposed. 

As another alternative, the trapezoidal' waveguides 106/106' could be replaced 
by the transition sections 130, shown in Figure 30, with or without spacings d } and d 2 , 
to form a periodic structure having sinusoidally- varying sides. It should be noted that 

10 these periodic structures are merely examples and not intended to provide an exhaustive 
detailing of all possibilities; various other periodic structures could be formed from 
unit cells comprised of all sorts of different shapes and sizes of elements. 

It should be noted that voltages can be applied to some or all of the strips in 
order to establish charges on the strips of the unit cells, which would change their 

15 permittivity and thus vary the optical transfer function of the periodic structures. If the 
dielectric material surrounding the strip is electro-optic, then the applied voltages 
would also change the permittivity of the dielectric, which also contributes to changing 
of the optical transfer function of the periodic structure. 

A two-dimensional photonic bandgap structure can be created by placing two or 

20 more of the periodic structures side-by-side to form a two-dimensional array of unit 
cells (comprised of strips of various shapes and sizes). A three-dimensional array can 
be created by stacking a plurality of such two-dimensional arrays in numerous planes 
separated by dielectric material. The size and shape of the strips are determined such 
that stop bands in the optical spectrum appear at desired spectral locations. 

25 Figure 40(a) illustrates an edge coupler 139 created by placing two strips 100" 

parallel to each other and in close proximity over a certain length. The separation S c 
between the strips 100" could be from about 1 |im (or less) to about 20 \xm and the 
coupling length L c could be in the range of a few microns to a few dozen millimeters 
depending on the separation S c , width and thickness of the strips 100", the materials 

30 used, the operating wavelength, and the level of coupling desired. Such a positioning 
of the strips 100" is termed "edge coupling". For example, if the strips are Au, the 
surrounding material is Si0 2 and the operating free space wavelength is near 1550 nm, 
good coupling lengths L c are in the range from about 100 to about 5 mm and good 
separations S c are in the range from about 2 |im to about 20 jim. Specific dimensions 

35 depend on the thickness and width of the strips and are selected to implement a desired 
coupling factor. 

The gaps between the input and output of the waveguide sections shown would 
ideally be set to zero and a lateral offset provided between sections where a change of 
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direction is involved. As shown in Figure 40(c), concatenated bend sections 108 
forming an S bend, as discussed in the embodiments associated with Figures 33 and 34 
could be used instead of the sections 104, 100 and 100" shown in Figure 40(a). Other 
observations made regarding the bend section disclosed in Figure 31 also hold in this 
5 case. 

Although only two strips 100" are shown in the coupled section, it should be 
understood that more than two strips can be coupled together to create an NxN. coupler. 

As illustrated in Figure 40(b) a voltage can be applied to the two edge-coupled 
sections 100" via minimally invasive electrical contacts. Figure 40(b) shows a voltage 

10 source 126 connected directly to the sections 100" but, if the sections 100, 104 and 
100" in each arm are connected together electrically, the source 126 could be 
connected to one of the other sections in the same arm. Applying a voltage in such a 
manner charges the arms of the coupler, which, according to the plasma model for the 
waveguide, changes its permittivity. If, in addition, the dielectric material placed 

15 between the two waveguides 100 is electro-optic, then a change in the background 
permittivity will also be effected as a result of the applied voltage. The first effect is 
sufficient to change the coupling characteristics of the structure but, if an electro-optic 
dielectric is also used, as suggested, then both effects will be present, allowing the 
coupling characteristics to be modified by applying a lower voltage. 

20 Figures 41(a) and 41(b) illustrate coupled waveguides similar to those shown in 

Figure 40(a) but placed on separate layers in a substrate having several layers 140/1, 
140/2 and 140/3. The strips could be placed one directly above the other with a thin 
region of dielectric of thickness d placed between them. Such positioning of the strips 
is termed "broadside coupling". The coupled guides can also be offset from broadside 

25 a distance S c , as shown in Figures 41(a)and 41(b). The strips could be separated by 
about d = 1 |im (or less) to about 20 ^m, the coupling length could be in the range of a 
few microns to a few dozens millimeters and the separation S c could be in the range of 
about -20 to about H-20 |im, depending on the width and thickness of the strips, the 
materials used and the level of coupling desired.- (The negative indicates an overlap 

30 condition) 

As before, curved sections could be used instead of the straight and angled 
sections shown in Figure 41(a). 

Gaps can be introduced longitudinally between the segments of strip if desired 
and a lateral offset between the straight and angled (or curved) sections could be 
35 introduced. 

Though only two strips are shown in the coupled section, it should be 
understood that a plurality of strips can be coupled together on a layer and/ or over 
many layers to create an NxN coupler. 
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As shown in Figure 41(b), a voltage source 126 could be connected directly or 
indirectly to the middle (coupled) sections 100" in a similar manner to that shown in 
Figure 40(b). 

It should be appreciated that in any of the couplers described with reference to 
5 Figures 40(a), 40(c) and 41(a), the adjusting means shown in Figure 40(b) could be 
replaced by either of the electrode configurations and biasing arrangements described 
with reference to Figures 36(b) and 36(c). It should also be noted that the electrode 
configuration and biasing arrangement could be applied to either or both of strips 100" 
of such couplers. 

10 As illustrated in Figure 42, an intersection 142 can be created by connecting 

together respective ends of four of the angled waveguide sections 104, their distal ends 
providing input and output ports for the device. When light is applied to one of the 
ports, a prescribed ratio of optical power emerges from the output ports at the opposite 
side of the intersection. The angles 4>i...4> 4 . can be set such that optical power input 

15 into one of the ports emerges from the port directly opposite, with negligible power 
transmitted out of the other ports. Any symmetry of the structure shown is 
appropriate. 

Various other modifications and substitutions are possible without departing 
from the scope of the present invention. For example, although the waveguide 

20 structure shown in Figures 1(a) and 1(b), and implicitly those shown in other Figures, 
have a single homogeneous dielectric surrounding a thin metal film, it would be 
possible to sandwich the metal film between two slabs of different dielectric material; 
or at the junction between four slabs of different dielectric material. Moreover, the 
multilayer dielectric material(s) illustrated in Figure 41(a) could be used for other 

25 devices too. Also, the thin metal film could be replaced by some other conductive 
material or a highly n- or p-doped semiconductor. It is also envisaged that the 
conductive film, whether metal or other material, could be multi-layered. 

Specific embodiments of modulation and switching devices 
30 Modulation and switching devices will now be described which make use of the 

fact that an asymmetry induced in optical waveguiding structures having as a guiding 

element a thin narrow metal film may inhibit propagation of the main long-ranging 

purely bound plasmon-polariton mode supported. 

The asymmetry in the structure can be in the dielectrics surrounding the metal 
35 film. In this case the permittivity, permeability or the dimensions of the dielectrics 

surrounding the strip can be different. A noteworthy case is where the dielectrics above 

and below the metal strip have different permittivities, in a manner similar to that 

shown in Figure 17(a). 
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A dielectric material exhibiting an electro-optic, magneto-optic, thermo-optic, 
or piezo-optic effect can be used as the surrounding dielectric medium. The modulation 
and switching devices make use of an external stimulus to induce or enhance the 
asymmetry in the dielectrics of the structure. 
5 Figures 43(a) and 43(b) depict an electro-optic modulator comprising two metal 

strips 110 and 112 surrounded by a dielectric 114 exhibiting an electro-optic effect. 
Such a dielectric has a permittivity that varies with the strength of an applied electric 
field. The effect can be first order in the electric field, in which case it is termed the 
Pockels effect, or second order in the electric field (Kerr effect), or higher order. 

10 Figure 43(a) shows the structure in cross-sectional view and Figure 43(b) shows the 
structure in top view. The lower metal strip 1 10 and the surrounding dielectric 1 14 
form the optical waveguide. The lower metal strip 1 10 is dimensioned such that a 
purely bound long-ranging plasmon-polariton wave is guided by the structure at the 
optical wavelength of interest. Since the "guiding" lower metal strip 110 comprises a 

15 metal, it is also used as an electrode and is connected to a voltage source 1 16 via a 
minimally invasive electrical contact 118 as shown. The second metal strip 112 is 
positioned above the lower metal strip 1 10 at a distance large enough that optical 
coupling between the strips is negligible. It is noted that the second strip also be placed 
below the waveguiding strip instead of above. The second strip acts as a second 

20 electrode. 

The intensity of the optical signal output from the waveguide can be varied or 
modulated by varying the intensity of the voltage V applied by the source 1 16. When 
no voltage is applied, the waveguiding structure is symmetrical and supports a 
plasmon-polariton wave. When a voltage is applied, an asymmetry in the waveguiding 

25 structure is induced via the electro-optic effect present in the dielectric 1 14, and the 
propagation of the plasmon-polariton wave is inhibited. The degree of asymmetry 
induced may be large enough to completely cut-off the main purely bound long-ranging 
mode, thus enabling a very high modulation depth to be achieved. By carefully laying 
out the electrodes as a microwave waveguide, a high frequency modulator (capable of 

30 modulation rates in excess of 10 G bit/s) can be achieved. 

Figures 44(a) and 44(b) show an alternative design for an electro-optic 
modulator which is similar to that shown in Figure 43(a) but comprises electrodes 
1 12A and 1 12B placed above and below, respectively, the metal film 1 10 of the optical 
waveguide at such a distance that optical coupling between the strips is negligible. 

35 Figure 44(a) shows the structure in cross-sectional view and Figure 44(b) shows the 
structure in top view. A first voltage source 1 16A connected to the metal film 1 10 and 
the upper electrode 112A applies a first voltage V! between them. A second voltage 
source 116B connected to metal film 110 and lower electrode 112B applies a voltage 
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V 2 between them. The voltages Vj and V 2 , which are variable, produce electric fields 
E x and Ej in portions 1 14A and 1 14B of the dielectric material. The dielectric material 
used exhibits an electro-optic effect. The waveguide structure shown in Figure 44(c) is 
similar in construction to that shown in Figure 44(a) but only one voltage source 1 16C 
5 is used. The positive terminal (+) of the voltage source 1 16C is shown connected to 
metal film 1 10 while its negative terminal (-) is shown connected to both the upper 
electrode 112A and the lower electrode 112B. With this configuration, the electric 
fields E { and produced in the dielectric portions 114A and 114B, respectively, are in 
opposite directions. Thus, whereas, in the waveguide structure of Figure 44(a), 

10 selecting appropriate values for the voltages V } and V 2 induces the desired asymmetry 
in the waveguide structure of Figure 44(c), the asymmetry is induced by the relative 
direction of the electric field above and below the waveguiding strip 1 10, since the - 
voltage V applied to the electrodes 1 12A and 1 12B produces electric fields acting in 
opposite directions in the portions 1 14A and 1 14B of the dielectric material. 

15 The structures shown in Figures 44(a), (b) and (c) can operate to very high 

frequencies since a microwave transmission line (a stripline) is in effect created by the 
three metals. 

Figure 45 shows in cross-sectional view yet another design for an electro-optic 
modulator. In this case, the metal film 110 is embedded in the middle of dielectric 

20 material 114 with first portion 1 14D above it and second portion 1 14E below it. 
Electrodes 1 12D and 1 12E are placed opposite lateral along opposite lateral edges, 
respectively, of the upper portion 114D of the dielectric 114 as shown and connected to 
voltage source 1 16E which applies a voltage between them to induce the desired 
asymmetry in the structure. Alternatively, the electrodes 1 12D, 1 12E could be placed V 

25 laterally along the bottom portion 1 14E of the dielectric 1 14, the distinct portions of 
the dielectric material still providing the asymmetry being above and below the strip. 

Figure 46 shows an example of a magneto-optic modulator wherein the 
waveguiding strip 1 10 and overlying electrode 1 12F are used to carry a current I in the 
opposite directions shown. The dielectric material surrounding the metal waveguide 

30 strip 1 10 exhibits a magneto-optic effect or is a ferrite. The magnetic fields generated 
by the current / add in the dielectric portion between the electrodes 1 10 and 1T2F and 
essentially cancel in the portions above the top electrode 1 12F and below the 
waveguide 110. The applied magnetic field thus induces the desired asymmetry in the 
waveguiding structure. The electrode 1 12F is placed far enough from the guiding strip 

35 110 that optical coupling between the strips is negligible. 

Figure 47 depicts a thermo-optic modulator wherein the waveguiding strip 110 
and the overlying electrode 1 12G are maintained at temperatures T 2 and T } 
respectively. The dielectric material 114 surrounding the metal waveguide exhibits a 
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thermo-optic effect. The temperature difference creates a thermal gradient in the 
dielectric portion 1 14G between the electrode 1 12G and the strip 1 10. The variation in 
the applied temperature thus induce the desired asymmetry in the waveguiding 
structure. The electrode 1 12G is placed far enough from the guiding strip 1 10 that 
5 optical coupling between the strips is negligible. 

It should be appreciated that the modulator devices described above with 
reference to Figures 43(a) to 47 may also serve as variable optical attenuators with the 
attenuation being controlled via the external stimulus, i.e. voltage, current, 
temperature, which varies the electromagnetic property. 

10 Figures 48, 49 and 50 depict optical switches that operate on the principle of 

"split and attenuate". In each case, the input optical signal is first split into N outputs 
using a power divider; a one-to-two power split being shown in Figures 48, 49 and 50. 
The undesired outputs are then "switched off or highly attenuated by inducing a large 
asymmetry in the corresponding output waveguides. The asymmetry must be large 

15 enough to completely cut-off the main purely bound long-ranging mode supported by 
the waveguides. The asymmetry is induced by means of overlying electrodes as in the 
waveguide structures of Figures 43, 46 or 47, respectively. 

In the switches shown in Figures 48, 49 and 50, the basic waveguide 
configuration is the same and comprises an input waveguide section 120 coupled to two 

20 parallel branch sections 122A and 122B by a wedge-shaped splitter 124. All four 
sections 120, 122A, 122B and 124 are co-planar and embedded in dielectric material 
126. The thickness of the metal film is d 3 . Two rectangular electrodes 128A and 
128B, each of thickness d l5 are disposed above branch sections 122A and 122B, 
respectively, and spaced from them by a thickness d 2 of the dielectric material 126 at a 

25 distance large enough that optical coupling between the strips is negligible. Each of 
the electrodes 128A and 128B is wider and shorter than the underlying metal film 122A 
or 122B, respectively. In the switch shown in Figure 48, the asymmetry is induced 
electro-optically by means of a first voltage source 130 A connected between metal film 
122A and electrode 128A for applying voltage V t therebetween, and a second voltage 

30 source BOB connected between metal film 122B and electrode 128B, for applying a 
second voltage V 2 therebetween. In the switch shown in Figure 49, the asymmetry is 
induced magneto-optically by a first current source 132 A connected between metal film 
122A and electrode 128A, which are connected together by connector 134A to 
complete the circuit, and a second current source 132B connected between metal film 

35 . 122B and electrode 128B, which are connected together by connector 134B to complete 
that circuit. 
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In the switch shown in Figure 50, the asymmetry is induced thermo-optically by 
maintaining the metal strips 122 A and 122B at temperature T 0 and the overlying 
electrodes 128A and 128B at temperatures T, and T 2 , respectively. 

It will be appreciated that, in the structures shown in Figures 48, 49 and 50, the 
5 dielectric surrounding the metal strip will be electro-optic, magneto-optic, or thermo- 
optic, or a magnetic material such as a ferrite, as appropriate. 

In general, any of the sources, whether of voltage, current or temperature, may 
be variable. 

Although the switches shown in Figures 48, 49 and 50 are 1x2 switches, the 
10 invention embraces lxN switches which can be created by adding more branch 
sections and associated electrodes, etc. 

It will be appreciated that, where the surrounding material is acousto-optic, the 
external stimulus used to induce or enhance the asymmetry could be determined by 
analogy. For example, a structure similar to that shown in Figure 45 could be used 
15 with the electro-optic material replaced by acousto-optic material and the electrodes 
1 12D and 1 12E used to apply compression or tension to the upper portion 1 14D. 

To facilitate description, the various devices embodying the invention have been 
shown and described as comprising several separate sections of the novel waveguide 
structure. While it would be feasible to construct devices in this way, in practice, the 
20 devices are likely to comprise continuous strips of metal or other high charge carrier 
density material, i.e. integral strip sections, fabricated on the same substrate. 

The foregoing examples are not meant to be an exhaustive listing of all that is 
possible but rather to demonstrate the breadth of application of the invention. The 
inventive concept can be applied to various other elements suitable for integrated optics 
25 devices. It is also envisaged that waveguide structures embodying the invention could 
be applied to multiplexers and demultiplexers. 

Although embodiments of the invention have been described and illustrated in 
detail, it is to be clearly understood that the same is by way of illustration and example 
only and not to be taken by way of limitation, the spirit and scope of the present 
30 invention being limited only by the appended claims. 
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